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Experiments and Impedance Modeling of Liners Including The Effect of Bias Flow 


Juan Fernando Betts 


(ABSTRACT) 


The study of normal impedance of perforated plate acoustic liners including the effect of 
bias flow was studied. Two impedance models were developed, by modeling the internal 
flows of perforate orifices as infinite tubes with the inclusion of end corrections to handle 
finite length effects. These models assumed incompressible and compressible flows, 
respectively, between the far field and the perforate orifice. 

The incompressible model was used to predict impedance results for perforated plates 
with percent open areas ranging from 5% to 15%. The predicted resistance results 
showed better agreement with experiments for the higher percent open area samples. The 
agreement also tended to deteriorate as bias flow was increased. 

For perforated plates with percent open areas ranging from 1% to 5%, the compressible 
model was used to predict impedance results. The model predictions were closer to the 
experimental resistance results for the 2% to 3% open area samples. The predictions 
tended to deteriorate as bias flow was increased. 

The reactance results were well predicted by the models for the higher percent open area, 
but deteriorated as the percent open area was lowered (5%) and bias flow was increased. 

A fit was done on the incompressible model to the experimental database. The fit was 
performed using an optimization routine that found the optimal set of multiplication 
coefficients to the non-dimensional groups that minimized the least squares slope error 
between predictions and experiments. The result of the fit indicated that terms not 
associated with bias flow required a greater degree of correction than the terms associated 
with the bias flow. This model improved agreement with experiments by nearly 15% for 
the low percent open area (5%) samples when compared to the unfitted model. The fitted 
model and the unfitted model performed equally well for the higher percent open area 
(10% and 15%). 
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List of Variables 


c speed of sound in the normal incidence tube 

c c speed of sound in the cavity 

Cd discharge coefficient (steady flow) 

Cdrag grazing flow drag coefficient 

d perforate hole diameter 

f frequency 

G nonlinear Bernoulli parameter 

H orifice inertial length parameter 

I complex amplitude of incident acoustic wave 

I.L. insertion loss 

J momentum of the bias flow jet 

Jo Bessel function of zeroth order 

J i Bessel function of first order 

J2 Bessel function of second order 

k wave number (=co/c) 

k r wave number in the r-direction 

L length of cavity 

M Mach number in the cavity 

M c bias flow Mach number in the cavity 
M g f mean grazing flow Mach number (=v gf /c) 

Mh bias flow Mach number in the holes of the perforate plates 

Mh effective Mach number in the holes of the perforate plates 

Msa acoustic Mach number in the holes of the perforate plates 

NIT normal incidence tube 

p pressure in the fluid 

^p 2 \ mean squared pressure 

p a acoustic pressure in the normal incidence tube 
Pah acoustic pressure in perforate’s orifice 
Pi incident acoustic wave on the liner 
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reflected acoustic wave from the liner 

far field pressure 

pressure in the perforate hole 

complex amplitude of reflected acoustic wave 

reflection coefficient 

reflection factor 

area of the boundary layer in perforate hole 

perforate hole area 

cavity area 

perforate thickness 

period 

fluid particle velocity in the x direction 
fluid particle velocity in the y direction 
fluid particle velocity 

acoustic particle velocity in the normal incidence tube 
acoustic particle velocity in perforate’s orifice 
bias flow velocity in the normal incidence tube 
bias flow velocity in the perforate’s orifice 
effective velocity in the normal incidence tube 
effective velocity in the perforate’s orifice 
effective velocity through orifice due to grazing flow 
mean grazing flow velocity 
rms acoustic particle velocity 
velocity in the boundary layer of the perforate hole 
fluid particle velocity in the z direction 


^W Iy ^ incident acoustic intensity to the liner 


(\V Ty ^ reflected acoustic intensity from the liner 

Z dimensional impedance 

a absorption coefficient 

y normalized reactance 
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5* boundary layer displacement thickness 

0 angle of incident acoustic wave on the liner with respect to the normal incidence 
T cavity propagation constant 

y Cp/c v = 1.4 for air 

k w shear stress from flow in the perforate duct. 

|i absolute viscosity of fluid 

g effective viscosity 

0 normalized resistance 

p density of fluid (air) in NIT 

p c density of fluid (air) in the cavity 
Ph density of fluid in the holes of the perforate samples 

o percent open area 

x time 

•o kinematic viscosity (=g/p) 

to circular frequency (=27if) 

E, normalized impedance 

cavity characteristic impedance 

Throughout this dissertation the word impedance means normalized impedance 
with respect to pc and the sign convention e +l<OT is used unless otherwise indicated. 
The words predicted and numerical are used interchangeably throughout this 
document. 



1. Introduction 


1.1 Objective 

Due to the current and projected concerns about community noise annoyance 
from commercial air traffic, a number of innovative noise reduction concepts are being 
considered. Acoustic liner treatments for engine nacelles have been an effective means 
of suppressing turbomachinery noise for over three decades. In the search for even 
greater liner efficiency, due in-part to the need to suppress turbomachinery noise 
emanating from high bypass-ratio engines, in-situ (in place) control of liner impedance 
has been an on-going goal starting with the work of Dean who explored bias flow as a 
means of changing liner impedances. Bias flow is the introduction of airflow, blowing or 
suction, perpendicular to the acoustic liner as seen in Fig. 1-1. Figure 1-2 shows the 
location of acoustic liners in a turbofan engine. 



Figure 1-1. Typical double degree bias flow liner configuration 1 . 
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PRIMARY ACOUSTIC ELEMENTS 



Bias flow can be viewed as a mechanism for optimizing a liner either for 
developmental purposes in a scale model test or perhaps eventually for an operational 
full-scale engine. 3 ' 4 Dean conducted an initial proof-of-concept study that showed several 
potential advantages 5 . The ability to control liner impedance “in-situ” (while the engine is 
operating) would allow several desirable possibilities including optimizing liner 
impedance to match different operating conditions or to more accurately match design 
conditions that could not be reached due to manufacturing tolerances. 

There are four objectives to this study: 

1. Acquire a quality experimental database of educed normal incidence 
impedances for perforates plates with and without bias flow 

2. Develop improved impedance models for perforated plates in the presence 
of bias flow 

3. Evaluate the models developed against the experimental database 

4. Use the experimental database to improve the model’s performance and 
evaluate the resultant model performance 

To this end the dissertation is divided into nine chapters and eight appendices. Chapter 2 
covers previous impedance modeling theory for perforated plates. Chapter 3 introduces 
new bias flow models, evaluation criteria, and other theoretical developments. This 
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chapter address objective 2. Chapter 4 describes the methods used to acquire the 
experimental database and therefore addresses objective 1. Chapter 5 shows impedance 
results of perforated plates without bias flow. This chapter is intended to analyze the 
performance of previous impedance models and to evaluate them against the evaluation 
error criteria to be used against the bias flow models. Chapters 6 and 7 present the 
impedance results of the experiments and two bias flow models developed assuming 
incompressible and compressible flow assumptions. This chapter addresses objective 3. 
Chapter 8 presents an experimentally fitted bias flow model. This chapter addresses 
objective 4. 

1.2 Approach 

An experimental database was produced which included perforate samples tested 
with and without bias flow in a single-degree-of-freedom (SDOF) liner configuration. 
The main trust of this dissertation work was the development and evaluation of bias flow 
impedance models for perforated plates. The database was used to evaluate the developed 
bias flow impedance models. After such evaluation was accomplished the models were 
“improved” through semi-empirical means to produce the best available frequency 
domain impedance model that includes the effect of bias flow. This approach permitted 
an assessment of how well a model “without” empiricism can predict experiments, and 
the relative improvement in the model’s performance with empirical corrections included 
in the model. 

The surface impedance predictions of the perforate-cavity system were obtained 
by way of a modified version of the NASA Langley Zwikker-Kosten Transmission Line 
Code (ZKTL) 6 . This computer program is based on Zwikker and Kosten’s theory for 
sound propagation in channels 7 . In general, the model is composed of continuous arrays 
of multi-degree-of-freedom liner elements. For the geometry being modeled, the 
perforated plate model is the uppermost element of the continuous lumped element array, 
as seen in Fig. 1-3. 

Matrix techniques are employed to compute the composite impedance due to the 
liner elements 8 . For the geometry of interest, these matrix techniques calculate the cavity 
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impedance (see Figure 1-3). Since flow is present in the cavity, wavenumber corrections 
were required and are presented in Section 3.5. Therefore, the total surface impedance is 
the sum of the cavity and perforated plate impedance. 



Perforated 

Plate 


Cavity 


Cavity 

Termination 

Boundary 

Condition 


Figure 1-3. Domain geometry being modeled and tested. 


In Chapter 3, new bias flow impedance models are developed for perforated 
plates. These impedance models are nonlinear, meaning that they are a function of the 
acoustic particle velocity. This can be expressed as 

S = — = f(v a ) (1-1) 

pcv a 

where p a , v a , p, c, and f(v a ) are the normalized impedance, acoustic pressure, acoustic 
particle velocity, fluid density, speed of sound, and an arbitrary function with respect to 
v a , respectively. Since the impedance in Eq. (1-1) is nonlinear, an iteration scheme was 
required, given the incident sound-pressure level (SPL), to make both sides of Eq. ( 1 - 1 ) 
equal. 
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In this study, the NASA Langley Normal Incidence Tube (NIT) was treated as a 
distributed element (channel) and the liner sheets with the backing cavity as elements. 
The different perforate impedance models were implemented in the lumped element part 
of the modular structure of ZKTL. The cavity backing plate was made porous to allow 
passage of bias flow, but with a high acoustic resistance to make it highly reflective 
(“hard”) as possible as seen in Fig. 1-3. 

Least squares linear fits (with zero intercept) between the experiments and 
predicted results were performed. The slope of the fit indicated the difference (systematic 
error) between experimental and predicted values, and the correlation indicated whether 
they were following the same trend. A slope and correlation of one indicated a perfect 
match between experiments and predictions. Acceptable limits on the deviations of slope 
and correlation from one will be used as a basis for establishing the evaluation criteria for 
the various models. 

Experiments were performed to determine the impedance of single degree of 
freedom liners with bias flow. The experimental data was obtained from the NIT. These 
liners were composed of nonlinear perforate facesheets followed by a 1.7 inch (4.32 cm) 
cavity and a high resistance fibermetal backing. Impedance data was acquired to 
determine the resistance of the fibermetal. It was determined that it had a minimal effect 
on the impedance of liner material. The bias flow was fed into a plenum chamber 3 
inches (7.62 cm) in length with a cross sectional area 2x2 inches (5.08x5.08 cm) before 
flowing through the high resistance fibermetal. 

Normal incidence impedance was measured using three stationary microphones. 
The first microphone was used to set the reference total sound pressure level (SPL) at the 
surface of the sample, and the other two microphones measured the transfer function 
between two points on the standing wave produced by the superposition of incident and 
reflected acoustic waves generated from the acoustic drivers and reflected from the 
perforate sample. The transfer function was then used to calculate the overall impedance 
of the sample-cavity system. 

The perforate samples were discovered to exhibit structural resonances in the 
middle to upper range of the frequencies tested. A vibration inhibiting post was inserted 
through the fibermetal to support the perforate sample in the center. Preliminary data 
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suggested that the post support eliminated the resonance below 3 kHz, which was the 
frequency range of interest. A thin nut was used to constrain the fibermetal that also had 
resonant frequencies within the range tested. Other than eliminating the structural 
resonance, the data showed that the post and the nut supports had a minimal effect on the 
impedance measurements. This determination was achieved by observing the continuity 
in resistance and reactance of the experimental impedance results. 

The perforate samples tested varied in open area from 1 - 15% with thickness to 
diameter ratios from 0.71 to 1.8. Twenty-three perforate samples were tested in all. The 
bias flow velocities tested ranged from 0 to 600 cm/s in the cavity. Tests were conducted 
over a frequency range from 1000 to 3000 Hz, one tone at a time (i.e., did not study 
frequency spectrum effects), in increments of 100 Hz. The reference SPL was set at 120, 
130, and 140 dB for low flow rates and at 130 dB at high flow rates where changing SPL 
had no effect on measured impedance (based on a few measurements over the SPL 
range). 

Acoustic resistance and reactance have been acquired for all samples. Overall, 
resistance increased with bias flow for all samples. At zero and low bias flow, increasing 
the reference SPL increased the resistance. Above a certain critical bias flow velocity; 
changing the reference SPL had no effect on the resistance. Therefore, the samples are 
said to exhibit linear behavior with respect to SPL at this critical bias flow rate. The 
reactance was minimally affected at low bias flow rates, but at high bias flow rates the 
reactance was significantly reduced. The effect was most noticeable as the velocity in the 
holes approached the choked condition. 
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2. Theory 


2.1. No Bias Flow Impedance Models 

2.1.1 GE Impedance Model 

Description and Major Assumptions 

• Incompressible Flow 

• Linear Term derived assuming Poiseuille Flow (Fully Developed DC Duct Flow) 

• Nonlinearity derived from non-dimensional analysis and empirical corrections 



Figure 2-1. The velocity profile in the hole of the perforate associated with the linear 
resistance term in the GE model. 

Motsinger and Kraft have discussed in detail the development of a Single-Degree- 
of-Freedom (SDOF) impedance model for perforate plates. 1 In their development they 
assumed a DC flow resistance equivalence with acoustic resistance, and therefore no 
frequency dependence in their impedance model, as seen in Fig. 2-1. Their resistance 
term was composed of three components; a linear viscous flow term, a nonlinear 
turbulent mixing loss component, and a grazing flow term. The GE resistance term, as 
described by Motsinger and Kraft, is 


0 


GE 


32t)t 1 

caC D d 2 2c(cC D f V> 


2 + 1.256— 
d 


r> 


M 


gf 


( 2 - 1 ) 


The linear resistance term in Eq. (2-1) can be derived by assuming Poiseuille flow 
in a duct, which for a Newtonian fluid creates a parabolic velocity profile within the duct 
as seen in Fig. 2-1. The velocity profile in the model is allowed to vary only in the radial 
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direction and the profile is fixed through the length of the duct. Motsinger and Kraft 
derived the nonlinear term through a non-dimensional analysis and an experimental fit. 
The GE mass reactance is given by 


— — — [t+ e d] (2-2) 

oC D l 


where e is an experimentally determined end correction 


0.85(1 -0.7Vq) 
l + 305(M gf J 


C D in the denominator of Eq. (2-2) was not included by Motsinger and Kraft, but it does 
appear when derived from a formal theoretical derivation. Kraft, Yu, and Kwan later 
corrected this problem. 2 Equation (2-2) can be derived by taking the one-dimensional 
momentum equation along the length of the duct and assuming the following: no 
viscosity, the velocity profile is constant through the cross section (i.e., “plug” flow), and 
the pressure field is only a function of the axial variation in the velocity field. 

Notice that the set of assumptions for the resistance and reactance are different. 
Some of these conflicts in assumptions as well as the omission of the discharge 
coefficient in the model may have been alleviated through the empiricism in the 
development of the other previously described model parameters. The complete GE 
model using Eqs. (2-1) through (2-3) becomes 


5 = 


32m 

ccC D d 2 


- + 


2c (cCj 


v, + ■ 


2 + 1.256 


V 


If 


M„ r +i- 


cC, 


-[t+e d] 


2.1.2 Crandall Impedance Model 

Description and Major A ssumptions 
• Incompressible Flow 


(2-4) 
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• Linear term derived assuming velocity profile has only radial variations and is 
constant through the length of the duct 

• Interaction effect between the holes included 

• Nonlinearity derived by assuming axial and radial variation in the velocity profile 

• The velocity is assumed to have a simple harmonic time dependence (i.e., there are no 
higher harmonics produced by the nonlinearity) 

• Acoustic parameters are replaced by their time averaged parameters for the nonlinear 
term 

Melting analyzed and discussed in detail the derivation of Crandall’s theory of 
acoustic propagation in perforates . 3 Crandall’s model assumes an infinitely long duct, and 
end corrections are added to account for the finite length of the perforate’s holes. 
Melting’s analysis is more of a theoretical nature, and relies less on empiricism than 
Motsinger and Kraft. His analysis yielded the following impedance model 



Here Jo and J| are the zero and first order Bessel functions. 


k' = 



( 2 - 7 ) 


is Stokes wave-number for a with highly conducting walls, and: 
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( 2 - 8 ) 



is Stokes wave-number for a wall with thermally nonconducting walls. The Fok function 
\| / is 



v|/(o) = Xa n (^)" 

n=0 


where 



ao=1.0 

ai=- 1.4092 

8 

ii 

o 

o 

a 3 =0.33818 

s 

II 

o 

o 

a 5 =0.06793 

a6=-0. 02287 

a 7 =0.003015 

a 8 =-0.01614 


The Fok function accounts for the acoustic interaction between the holes. According to 
Melling, Fok derived this solution of the impedance coupling between the perforate s 
holes for an infinitely thin plate. 

Besides the acoustic interaction of the holes, the Crandall impedance model 
introduces frequency dependence to the acoustic impedance. The linear frequency 
dependence was lacking in the analysis by Motsinger and Kraft. The frequency 
dependence in the linear term of the Crandall model is due to the Bessel function 

solution. 

Historically, Equation (2-5) was not normally used to calculate the impedance of 
a perforate, due to the difficulty in calculating the Bessel function with a complex-valued 
argument and separating the real and imaginary components of the impedance function. 
Although this impediment has been a problem in the past, current computers can easily 
handle the above-mentioned problems. However, the low and high frequency 
approximations to Eq. (2-5) are widely used, to provide additional physical insights to the 

Crandall impedance function. 

The low frequency approximation (also called the Poiseuille model) is 
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( 2 - 10 ) 


d 

2 



< 


1 


Making this approximation, Eq. (2-5) becomes 


32t)t | (l-o 2 ) . k (A 8d " 

coC D d 2 2c(aC D ) 2 3 aC D ^3 37 c\|/'(ct) ^ 


( 2 - 11 ) 


The first term in Eq. (2-11) is the same as the viscous DC linear term presented by 
Motsinger and Kraft for the GE model. Note that the first term in the reactance is 1/3 
larger in the Crandall model than in the GE model. The effective mass in the Crandall 
model is larger because the viscosity within the hole increases the effective mass within 
the hole. This is due to the fact that the velocity profile is parabolic in the Crandall 
model, and not “plug-flow” as assumed for the reactance derivation in the GE model. The 
lower average velocity in the parabolic distribution leads to a higher effective mass 
because of conservation of momentum. 

The high frequency approximation (also called the Helmholtz model) is valid if 


d 

2 



>10 


( 2 - 12 ) 


Making this approximation Eq. (2-5) becomes 


, 2 . 82 ^- 14,1 

coC d d 2 c(oC d ) ^ocjCp coC d d 3n\)/(a) 


(2-13) 


This approximation introduces frequency dependence in the linear term in the resistance 
that is absent in the GE model. Furthermore, the frequency independent term of the 
Poiseuille model vanishes. Also note that the effective mass is no longer a constant 1/3 
higher, but contains an additional attached mass that is frequency dependent. 


12 



In practice neither the low nor high frequency approximations are used, because 
of the desire for a single model that works for all frequencies. Therefore, the high 
frequency model is usually corrected to account for low and intermediate frequencies. 
Kraft, Yu, and Kwan provide one such correction by adding half of the Poiseuille viscous 
term to the Helmholtz model. 2 This model is called the GE-Rohr model, and is given by 


4 = 


16ut 

caC D d 2 


Vcop t 

• + 2.82 + 


coC d d 


(!-<+) 

2c(aC,f 


v. +1 


J^ +2 . 82 i™i + .l * 

oC D coC D d 3 7t V(a) 


A 


(2-14) 


2.1.3 Hersh 75 Impedance Model 

Description and Major Assumptions 

• In the near field the flow moves in the radial direction and is axisymmetric. 

• Within the near field the flow is incompressible and unsteady. 

• At low sound pressure levels the linear regime dominates 

• At high sound pressure levels the nonlinear regime dominates 

• Core flow dominated by entrance effects 


Displacement 

Boundary 

Layer 



Figure 2-2. Velocity profile for Hersh 75 model for a perforate hole. 

In 1975 Hersh and Rogers followed a different approach to modeling the 
impedance of the perforate from an impinging acoustic field. 4 According to them, 
perforate plates are too thin to assume fully developed duct flows through them as both 
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the linear components of the GE and Crandall models do. Hersh and Rogers concluded 
that since the orifice diameters are small in perforates, end corrections are not 
corrections, but a major component of the model’s viscous losses. Major viscous losses 
do not occur inside the perforate’s hole, but instead occur at the surface of the perforate. 
Consequently, their analysis focuses on modeling the near field of the perforate’s hole 
rather than the interior, as seen in Fig. 2-2. 

Hersh and Rogers developed their model by considering two regions, which they 
called the linear regime and nonlinear regime for low and high sound pressure levels 
respectively. In their model the linear and nonlinear regimes are determined by 


acod 


(2-15) 


If this ratio is less or greater than one, the regime is linear or nonlinear respectively. 
For the linear regime the impedance is 


where 


^ Hersh 75 


= k(d + 1) 


qcod ! L 


3tcC 


ifi+4l ^ 


1 + J— - 

Re 


8 


acod 


(2-16) 
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Nonlinear loss Viscous loss 
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V2Re 
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Viscous loss 


* S K)1 
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Nonlinear loss 


(2-17) 


Re = 


co(d + 1) 2 


D 


(2-18) 
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The resistance term has two dominant terms, a nonlinear and a viscous loss term. These 
terms appear in both the resistance and reactance, indicating a linear and nonlinear 
coupling between resistance and reactance. This coupling is dominated by the viscous 
term when Ratio (2-15) is less than one. Rewriting Eq. (2-16) and (2-17) with Ratio (2- 
15) much less than one produces 


0 


Hersh75 


k(d + t) 



1 

H 

Re 


\ 

i 

+ ik(d + t) 
y 

+ sl2Re_ 


(2-19) 


Since the impedance is dominated by the viscous term, it is directly related to the 
displacement boundary layer along the surface of the perforate. Furthermore, the effects 
of viscosity increase the mass reactance. This result is in line with the Crandall model, 


where viscosity inside of the hole also increased the mass reactance. Therefore, the effect 
of viscosity increases the mass reactance both in the near field and within the hole. 

For the nonlinear regime, Ratio (2-15) is larger than one. The impedance becomes 


'Hersh75 


= k(d + 1) 



1 .03C n toq(d + 1) 

v. 


( 2 - 20 ) 


XHersh7J — ^(d + 0 


0.33 

c D 





( 2 - 21 ) 


In this regime nonlinear effects dominate the impedance, and the viscosity effects 
disappear in both the resistance and reactance. 
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2.1.4 Hersh 99 Impedance Model 

Description and Major Assumptions 

• Incompressible Flow 

• Inviscid core flow with a viscous boundary layer inside the perforate hole 

• Developing viscous boundary layer 

• Unknown assumed model parameters determined through experiment 

• t/d»l 



Figure 2-3. Velocity profile in the perforate for Hersh 99. 

Hersh, Walker, and Celano 5 developed another impedance model in 1999, where they 
assumed an inviscid core with a boundary layer profile within the hole, as seen in Fig. 2- 
3. They developed this model from a control volume approach using the principles of 
conservation of mass and momentum. An assumed velocity and viscous loss was then 
inserted into these integral equations. The resultant equations become a function of 
unknown parameters that are themselves a function of geometrical variables. Then the 
functional relationship between the unknown parameters in the integral equations and 
geometrical variables is determined through experiments. A more detailed analysis will 
be presented in Appendix A, where a bias flow model using Hersh’s ideas is presented. 
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2.2 Bias Flow Impedance Models 

2.2.1 Dean’s Bias Flow Model 

Description and Major Assumptions 

• Same assumptions as the Hersh 75 model 

• Acoustic particle velocity replaced by bias flow velocity 

• Additional impedance corrections determined empirically 

Bias flow has been a concept for “in-situ” control of impedance in acoustic liners. 
Dean conducted an initial proof of concept in 1976 6 and developed a bias flow model 
based on the principles outlined by Hersh in his 1975 model. Dean, like Hersh, divided 
the flow region into two regions, the linear and nonlinear regimes respectively. He 
utilized the same terms as Hersh (Ratio (2-15)) to determine which impedance region was 
operating (linear or nonlinear), except that he replaced the acoustic particle velocity with 
the bias flow velocity in that term. 

Therefore, the term that Dean used to determine the operating impedance regime 

was 

_A_ (2-22) 

acod 


He also replaced the acoustic particle velocity with the bias flow velocity in both the 
linear and nonlinear impedance model regimes. Therefore his model ignores the effect of 
sound pressure change caused by changes in the acoustic particle velocity. He assumes 
the bias flow velocity to be much greater than the acoustic particle velocity, and therefore 
SPL effects are considered negligible. If the ratio (2-35) is less than one, his impedance 
model is 
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(2-24) 


The last two terms in the resistance represent bias flow and grazing flow, respectively. 
Dean does not explain where these terms come from except to say “the effects of bias 
flow... has been experimentally investigated sporadically... and in broad terms the d-c 
flow resistance increases in direct proportion to the bias flow velocity and inversely to the 
square of the porosity.” It appears that Dean might be alluding to a similar empirically 
determined nonlinear effect described by Motsinger and Kraft in the GE model 
previously discussed in this chapter. Dean’s model has further differences to Hersh. 
Although the models are similar, they are not the same after accounting for the change in 
v a by Vb. Dean does not explain these differences. 

The nonlinear regime occurs when ratio (2-35) is greater than one as with Hersh. 
The impedance becomes 
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(2-25) 


(2-26) 


This impedance model is similar to the Hersh 75 model, and has the same bias and 
grazing flow corrections. Dean does not shed any more light into the development of this 
nonlinear regime impedance model. 


2.2.2 Previous Bias Flow Corrections 

The acoustic impedance models presented, except for Dean’s Bias Flow model, 
do not account for the effect of bias flow. As seen in Dean’s model, bias flow was 
accounted by replacing the acoustic particle velocity with the bias flow velocity and 
stating that the bias flow velocity was much greater than the acoustic particle velocity. As 
will be shown later, replacing the particle acoustic velocity with the bias flow velocity is 
questionable. 
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In the NASA Langley Zwikker-Kosten Transmission Line Code (ZKTL) 7 , bias 
flow is accounted for by replacing the acoustic particle velocity in the impedance models 
with 


v a + v b (2-27) 

This correction seems more reasonable than Dean’s correction, especially at low bias 
flow velocities. Nevertheless, this bias flow correction as well as Dean’s correction does 
not follow from basic principles. It is important to note that this correction has never been 
published, assessed or validated, and therefore was just a placeholder in ZKTL until a 
better correction was developed. 

In another study, Premo derived a bias flow model using a time-domain approach. 
This approach yielded the following replacement for the acoustic particle velocity to 
include bias flow 8 


V( 1 - 15v ™ s ) 2 + (2v b ) 2 (2-28) 

The relevance and comparison of this bias flow correction is further discussed in Section 
3.1. 


2.3 Grazing Flow 

The effect of grazing flow on the acoustic impedance of perforates has been 
extensively studied for the past three decades. Rice developed an empirical acoustic 
resistance model for a single hole is given by 9 

0 = 0.3M gf (2-29) 

Equation (2-42) was modified in the same reference to account for percent open area (ct) 
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( 3. 15x1 O' 2 


(2-30) 


+ 1.770 + 0.3 

Since the development of this empirical model, there have been several grazing flow 
impedance models developed. Rice devised a model that assumed that vortices were 
formed when the fluid interacted with the acoustic jet exiting the perforate hole . He 
presented an approximate simplified solution, where 


0oc 


dM gf 

08 


(2-31) 


This semi-empirical model showed that the grazing flow impedance was a function of 
boundary layer thickness (5). Boundary layer thickness was ignored in the empirical 
model in Ref. 9. Armstrong, Beckemeyer, and Olsen devised a method of relating the 
acoustic propagation constant k z along the length of the duct to the acoustic impedance in 
the presence of grazing flow using a “waveguide method. This method took into 
account the boundary layer thickness through the velocity profile of v g f. 

Rice developed a theoretical model for grazing flow, 12 using the principles 
employed in the impedance model outlined by Hersh and Rogers (Hersh 75 model, see 
Section 2.1.3) 4 . Rice solves a simplified form of the Navier-Stokes equations assuming 
only a radial component of velocity entering the hole in a spherical coordinate system, 
and included grazing flow in his analysis. The resultant grazing flow part of the model 

was 




(2-32) 


Heidelberg, Rice and Homyak modified Eq. (2-45) to account for the effect of the grazing 
flow boundary layer thickness (5) 13 . Their model is given as 
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(2-33) 


0 = 
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3. New Theoretical Developments 

There has been two major approaches to acoustic impedance modeling of 
perforate plates presented in Chapter 2. In the first approach, perforated plates were 
considered to be too thin for the orifices to be modeled as cylindrical ducts. Hersh and 
Sivian use this approach, where the end effects and boundary layer are significant. The 
second approach has been to model the perforate orifices as cylindrical ducts, an then to 
correct for “end effects.” Crandall and Kraft have implemented this approach. 

Although both of these approaches have been used to develop acoustic impedance 
models for perforates, the second approach has by far been more popular than the first. 
To date, the Crandall model as presented by Melling is the most comprehensive and 
complete impedance model study for perforated plates that 1 have found. This model 
contains fewer approximations, relies less on empiricism, and includes more effects than 
any of the other model studies. For this reason, the bias flow impedance models 
developed in this chapter will model the perforate orifices as cylindrical ducts. Chapters 6 
through 8 evaluate these models using experimental results. Appendix A shows the 
development of a bias flow impedance model using the approach outlined by Hersh. This 
model is not evaluated in this dissertation. 


3.1 Perforate Bias Flow (PBF) Model 

The momentum equation for a viscous fluid is 


P 


— + (v*V)v 

3t I 


-Vp + |iV 2 v 


(3-1) 


where the compressibility term J^|iV(V» v) is assumed small and therefore is omitted 

from Eq. (3-1). Note, that by making this approximation, density disturbances are 
ignored. Let any fluid variable q be the sum of a mean flow qb and an acoustic component 
q a , that is 
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q(x,r,x)= q b (x,r,x)+q a (x,r,x) 


(3-2) 


Substituting Eq. (3-2) into (3-1) produces 


~ dv a 

3t 


+ ((v b + v a )*VXv b + v a 


= "Vp b - Vp a +llV 2 v b + pV 2 v a 


(3-3) 


where 


v b = v v (x,r)e x 
v a = v a (x,r,t;)e x 


(3-4) 


and simplifying yields 


^ + (v b • V )v b + (v b • V)v a + (v a • V)v b + (v a • V)v a 
= "V Pb - V Pa + |iV 2 v b + pV 2 v a 


(3-5) 


This equation contains bias flow only, coupled bias flow/acoustic, and acoustic only 
components. The bias flow only components should balance, leaving only coupled and 
acoustic only components. Furthermore, assume that perforate holes can be modeled as 
cylindrical ducts, and that only a velocity component in the x-direction along the length 
of the perforate hole exists. This velocity component is assumed to vary in the x and 
radial directions. So far, the above description pertains to the interior of the perforate 
hole. To indicate this, the variables v a and v b will be subscripted with v a h and v bb . Making 
these approximations to Eq. (3-5) and letting v a h and v b h now represent the acoustic and 
bias flow velocity in the x direction inside the perforate’s hole produces 
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Rewriting Eq. (3-6) and grouping linear and nonlinear coupling terms yields 
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The right hand side of Eq. (3-7) can be written as 


^Pah __ 

r 3p ah 

+ dp ah 

\ 

3x 

[dx 

linear 9x 

nonlinear j 


Therefore the linear and nonlinear terms of Equation (3-7) balance their corresponding 
counterparts respectively, in Eq. (3-8). By utilizing the definition of impedance, the linear 
term of Eq. (3-8) can be rewritten in terms of the linear impedance by noting that 


dp 


ah 


dx 


linear 




(3-9) 


Substituting the above relationship in Eq. ( 3 - 7 ), assuming dx-t, dp a h=Pah, pah~Pa, and 


dv ah =v a h produces 

-^linearP CV »h-p( V bh V ah)-^ L = Pa ( 3 ‘ 10 ) 


Figure 3-1 depicts a control volume where the incident and hole variables are defined. 
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Perforate Hole 
Vah, V bh , S s 


Figure 3-1. Control volume for the incident and hole regions. 


From continuity, the relationship between the incident and hole velocities for an 
incompressible fluid is given by 


V a ~ ^D^ah 
V b = ^D^bh 


(3-11) 


Inserting the relationships in Eq. (3-11) into Eq. (3-10) and multiplying by v a produces 


- ^l' n «rP CV a P y y 2 P^a _ _ 

C D c b ■ 2(C 0 af P - 1 


(3-12) 


The following procedure was used by Melling to relate the acoustic particle velocity to 
the rms velocity. Integrating both sides of Eq. (3-12) over the period yields 


B l ,nearP CV a dT + lf 
0 C D a T J 0 (C D oJ 




T 


T 


JPa V adX 


(3-13) 


where v^ has been replaced by |v a j vf in the third term to ensure that it remains positive. 

This is done for the following reason. A harmonic solution of the form Cos(ot) is 
inserted for v a into Eq. (3-13) in every term except the third term. Rewriting the resulting 
equation with the third term on the right hand side yields 
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The left hand side of Eq. (3-14) is a constant regardless of the sign of v„ yet the right 
hand side is either positive or negative depending on the sign of v a (remember that a 
harmonic solution is assumed for v a ). This can be expressed mathematically as 


A = C 
and 

A = -C 


(3-15) 


where A and C are constants. The solution to Eq. (3-15) can be A=C=0, which is the 
trivial solution. Another solution is for the right hand side to be in absolute values. 
Therefore, in order to eliminate the trivial solution and for the equality to hold regardless 
of the sign of v„ v 3 a must be replaced by |v a |vj . Making this replacement into Eq. (3-13) 

produces 


5 lmearP£ v 2 + P 


C D o 




(3-16) 


Dividing both sides by pcvj^ produces 


1 -^ ssl + 1 -.[2v b + 1.2v„J 

C D o 2c(C D of 


(3-17) 


Melling notes that the pressure drop across a “sharp edge” orifice has been studied in 
some detail 1 . From these studies, a departure from the nonlinear term in Eq. (3-17) is 
suggested. Making the corrections indicated by Melling to Eq. (3-17) yields 
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[2v b + 1.2vJ 


(3-18) 
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The linear impedance (linear) term is determined by solving Eq. (3-7) utilizing only the 
linear terms. Writing that equation leaving only the linear term yields 
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The linear impedance of the Crandall model (see Section 2.1.2) is the harmonic solution 
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of Eq. (3-19) assuming that the term 
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8x 2 


0 . This simplification effectively assumes 


there are no acoustic waves traveling along the length of the duct, or the wavelength is 

9x 2 

harmonic solution of the following form is assumed: 


much greater than the perforate thickness since -^-y « 0(k 2 ) « 1 . To solve Eq. (3-19) a 


v ah (r,t) = v(r)e' (<0, > 
p ah (x,t) = p(x)e i(<M) 


(3-20) 


Inserting the above relations into Eq. (3-19) yields: 


d 2 v 1 dv 

r-H (■ 

dr 2 r dr 



(3-21) 
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Equation (3-21) is the inhomogeneous Bessel equation of zero order. Assuming the 
pressure gradient in Eq. (3-23) to be constant, the solution of this equation is 


v(r) = - 




J 0 [k s r] 

J 0 [k s R] 


(3-24) 


where the no-slip condition at the tube wall, v(R) = 0, has been imposed. R is the radius 
of the perforate hole. Averaging vover the cross-sectional area produces the mean 
velocity 


{*) ~ ” 7 Jv (r) ( 2 itr)dr 

o (3 .25) 
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From Eq. (3-23) Ap is: 


Ap = <{>t = 
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where 



(3-27) 


Using the fact that the acoustic impedance is 



and the diameter is twice the 


radius produces: 
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(3-28) 



Sivian showed that the viscosity inside the perforate hole along a highly thermally 
conductive wall was different from the absolute viscosity outside the hole 2 . Sivian 
determined that for air over a wide range of temperatures the relationship between jj. and 
fj.' is: 

|i=2.179p (3-29) 


Noting that the analysis for the impedance was done within the perforate’s hole k s 
becomes k' . Equation (3-28) is the first term of the Crandall model as shown in Section 

2.1.2. Tijdeman has studied extensively the propagation of sound waves in cylindrical 
tubes 3 . He presents results of acoustic propagation in ducts under less restrictive 
assumptions than what has been presented here. His analysis was not used for this study 
because a simpler model would provide a good initial understanding of the dynamics of 
the bias flow perforated plate. Furthermore, it was not known whether the approach of 
modeling the perforate orifices as ducts would work, and therefore given the resource 
constraints of this study a more complicated model was not warranted. Nevertheless, 
Tijdeman analysis could be used in future work as a starting point for developing a more 
comprehensive impedance model. 

Melling also describes contributions by Sivian and Ingard for end effects to the 
perforate hole, and Fok’s contribution to interaction effects between holes. His analysis is 
not going to be repeated here, but his results will be included. Therefore the Perforate 
Bias Flow model is: 



(3-30) 
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where the Fok function \ ;'(o) is defined in Eq. (2-9). 

Note that the Perforate Bias Flow model is similar to the Crandall Impedance 
model. The Crandall Impedance model can be “corrected” to become the Perforate Bias 
Flow model by letting 


V, =2v b + 1.2v nM (3-31) 

A problem arises when large enough negative bias flow rates would lead to 
negative impedance values. A recent study has experimentally shown the effect of 
negative bias flow rates of two degree-of-freedom liners. 4 The results of this suggest that 
increasing negative bias flow rates do not necessarily lead to lower or negative 
impedances. Consequently, Eq. (3-31) needs to be restricted to positive values. A 
possible method to accomplish this is to square Eq. (3-31) and take its square root. Doing 
this produces 


v a = A /(l. 2 v ms ) 2 + (2v b ) 2 +4.8v nns v b (3-32) 


Thus the sign of v b does affect the magnitude of impedance but cannot drive it negative. 
Therefore, blowing produces a different impedance result than suction. Comparing Eq. 
(3-32) to Premo’s bias flow correction (Eq. (2-41)), it is noted that Premo s correction 
lacks the acoustic-bias flow coupling term 4.8v rms v b and therefore that model predicts the 
same resistance whether bias flow is provided via blowing or suction. 

Equation (3-30) has a low and high frequency approximation. These 


approximations are valid if 



(3-33) 


respectively. The low and high frequency approximations of Eq. (3-30) are 
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respectively. The PBF model runs into the problem that the Bessel functions in the linear 
component of the impedance obscure its physical relationship to geometrical parameters 
such as o, t, and d. The low and high frequency approximations to the PBF model do not 
account for intermediate frequencies. Therefore, a single model that does not contain the 
Bessel functions, and “works” for all frequencies is desirable. One such correction similar 
to one done by Kraft, Yu, and Kwan, is 5 
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Equation (3-36) will be called the Perforate Bias Flow Intermediate Frequency model 
(PBFIF) model. This is the model that will be evaluated against experiments. The exact 
model presented in Eq. (3-30) would be more accurate than Eq. (3-36). Nevertheless, 
during the evaluation process, using the PBFIF model provides better insight into the 
reasons for the good or bad agreement between the model and experiments. Using the 
PBFIF model also provides an extra degree of freedom for the experimental fit of this 
model presented in Section 3.6. 


3.3 Perforate Compressible Bias Flow (PCBF) Model 

The models studied so far, had a key assumption. These models assumed 
incompressible flow through the perforate holes. When bias flow is added, Mach 
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numbers within the perforate plate’s hole can be substantial, and therefore fluid 

compressibility must be taken into account. 

Compressibility corrections can be added to the models by allowing the incident 
field mean density to be different from that of the perforate’s hole. During the derivation 
of both the PBF model and the PBLBF model the use of a control volume approach was 
employed. The control volume was used to relate the incident field acoustic and bias flow 
velocities to that of the perforate’s hole, respectively. From the continuity principle, the 
incident and hole flow values can be related to each other. Therefore, Eq. (3-11) becomes 


v a =— aC D v ah 

P (3-37) 

v b =— oC D v bh 
P 

where the mean hole density, p b needs to be determined. 

To determine the incident and mean perforate hole density, isentropic conditions 
can be employed. For an ideal gas assuming constant specific heat c p , the first law of 
thermodynamics can be written as 6 


cT = c T + — O 38 ) 

p o p 2 

where T c is the stagnation temperature. The speed of sound is given by 

c 2 = yRT C 3 ' 39 ) 

Inserting Eq. (3-39) into Eq. (3-38) and rewriting in terms of the incident and hole 
properties produces 
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For an isentropic fluid the following equation for a perfect gas can be written as 7 
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Inserting Eq. (3-41) into Eq. (3-40) and rearranging yields 
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The principle of continuity can be written as 

A = c Pi- v bh = C Ph M h c h 
S h ° P v b ° p Me 


(3-43) 


Substituting Eqs. (3-42) and (3-40) into Eq. (3-43) and assuming a circular jet produces 


o 


l + -(y-l)M 2 
1 + 1(Y-1)M 2 
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(3-44) 


Mh is computed through iteration, and p h is calculated from Eq. (3-42). 

The compressible momentum equation with the simplifications given by Eqs. (3- 
1) through (3-9) yields 
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(3-45) 
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Inserting Eq. (3-37) into the above relation, multiplying by v a and integrating over the 
period, noting that |v t | needs to be used in the third term to prevent it from being 

negative, produces 
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Following a similar derivation as the one outlined in Section 3.1 produces 
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where 
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Notice that the compressibility term y only appears in the nonlinear term of Eq. (3- 

/ rh 

47). The corresponding low and high frequency approximations are 
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respectively. An intermediate frequency model such as Eq. (3-36) is 
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(3-52) 


Equation (3-52) will be called the Perforate Compressible Bias Flow Intermediate 
Frequency (PCBFIF) model. 

3.4 Grazing Flow Model 


The following derivation uses some of the principles outlined by Rice 8 . It is 
important to note that although the general assumptions in this development follow from 
Rice, this model is more general and includes many more effects ignored by Rice. 
Moreover, this model includes bias flow, which was not included by Rice. Figure 3-1 
shows the geometry of the model for the grazing flow-bias flow perforate hole 
interaction. 
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Figure 3-1. Grazing flow bias flow perforate hole interaction. 

The jet emanating due to the bias flow interacts with the grazing flow. In the model that 
follows, the jet is treated as a rigid cone in the presence of grazing flow. The grazing flow 
generates vortices as it hits the cone. These vortices shed at a frequency given by the 
Strouhal number (St) 9 : 


<ad(y) 

” v g (y) 


( 3 - 53 ) 


Here d(y) and V g (y) are the diameter of the inviscid core, and the grazing flow velocity, 
respectively at some distance along the y-axis. V g (y) can be approximated by 

V(y) = ^I ( 3 - 54 ) 

where V g is the free field grazing flow velocity and 8 is the boundary layer thickness. The 
acoustic pressure generated by the vortices at point A is given by 


Pa = 


CpragPVg (WV™ 

2 


( 3 - 55 ) 
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where W is the distance that the jet emanates from the perforate hole into the grazing 
flow stream. To determine the radiated acoustic field, the convective wave equation is 
used. The convective wave equation in spherical coordinates assuming only radial 
variations is 


'3 ! p + 2 3p' 

II 

K r*J 



d 2 P „ 32 

— T + 2v r 

at 2 r a r at 


2 _\ 


a 2 P 2 a 2 P 

~ +v 'aT 2 


(3-56) 


The assumed solution to Eq. (3-56) is 


P = 


— e i ^“ ,-k ' r ) 
r 


(3-57) 


To determine the dispersion relationship, Eq. (3-57) is inserted into Eq. (3-56). This 
produces 


(M 2 -l) k2 (i + 2kM r ) h + |V 

r r r r 

V 


2M 


+ i- 


2k M 


= 0 


(3-58) 


and M r is the radial component of Mach number due to the bias flow. The roots of Eq. (3- 
58) are solved with respect to k to find the wave number corrections. Since the bias flow 
jet comes out with Mach number (M y ) in the y direction only, the radial component of 
this velocity needs to be computed to get M r . It can be shown from Fig. 3-2 that M r is 
related to M y by 


M =■ 


M v Vx 2 + W 2 


W 


(3-59) 


where M y is the bias flow Mach number leaving the perforate hole. M y is given by 


38 


(3-60) 


M 


y 


cc 



Figure 3-2. Geometry of Perforate Hole in Relation to the Acoustic Source. 
Therefore, the radiated pressure field is 


i( r ) = 


C Drae r o pVg(W) . r 


2r 


(3-61) 


where e iw has been omitted from Eq. (3-61) for clarity. At the perforate’s hole surface, r 
is given by 


r = ( 3 ' 62 ) 


Therefore the pressure field at the perforate hole surface is 

/ n C Drag r o pV E (W) ^_ ik j w i +t 2 (3-63) 

Pa X 2-n/w 2 +x 2 

That distance W is due to the sum of the acoustic and bias flow fluid particle 
displacements and is given by 
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(3-64) 


W = — + W t 
i(0 

where W b is the jet displacement into the grazing flow due to the bias flow. To determine 
W b the following analysis is made. The momentum and continuity equation for a circular 
jet can be written as 10 


y dv + 9u _ (i 1 9 f 3v 
dy dx p x dx l 3x 


(3-65) 


+ 3u + u 
9y 3x x 


(3-66) 


respectively. Velocity components along the x and y-axis are denoted by u and v, 
respectively. The y and x-axis are the lengthwise and radial components of the jet, 
respectively. The boundary conditions are 

3v 

x = 0 : u = 0; — = 0 

dx (3-67) 

X = oo ; v = 0 


The solution to Eqs. (3-65), (3-66), and (3-67) is 11 


3 K 1 
8rt uy (1 + 0.25P 2 ) 2 


(3-68) 


u = 0.25 



VfC p-0.25p 2 
y (l + 0.25p 2 J 


(3-69) 


where 
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(3-70) 



(3-71) 


J is the momentum of the jet given by 


J = 27 ipjv 2 xdx 

0 


(3-72) 


For the case of the perforate orifice, the equation for J becomes 



(3-73) 


At x=0, p=0, and v is 


0.09375d 2 v 2 

4uo 2 y 


(3-74) 


The distance the jet extends into the grazing flow stream is subjective. The fluid velocity 
continually slows down due to the viscosity in the surrounding fluid, and therefore where 
the jet is no longer affecting the grazing flow is debatable. Consequently, an assumption 
is made that the jet affects the grazing flow stream as long as it retains 10% of its original 
velocity. Inserting this assumption into Eq. (3-74) and simplifying yields 

y = W b = 0.2344^ (3-75) 
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Equation (3-63) must be space averaged to get an approximate impedance at the 
hole surface. Therefore the average acoustic pressure at the surface of the perforate due to 
grazing flow p g is 


Y 

P g =- 7 /p A ( x )dx (3-76) 

-Y2 

The total acoustic pressure at the surface of the perforate p tot is given by 

P,o, =P S +P g (3-77) 

where p s is another acoustic pressure source term present in the system. The bias flow 
impedance models presented previously give the acoustic particle velocity response. 
Therefore 


— = ^p+Sc.v (3-78) 

pcv a 

where £, p and ^ 3 v are the perforate bias flow impedance model and cavity impedance 
respectively. To solve Eq. (3-78) the left and right hand side of this equation must be 
equal. Therefore, through an iterative method an appropriate v a is found that makes both 
sides of this equation equal. 


3.5 ZKTL and Transmission Matrix Theory 

Jones and Parrott first presented the Zwikker-Kosten Transmission Line code 
(ZKTL) in 1995. 11 The idea of this code was to use transmission matrices to calculate 
pressure and velocity at each section for a multi-element liner. The forward transmission 
matrix through a dissipative channel has the following form 12 
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(3-79) 


P m+l 


T„ 

fO 

1 

1 

3 

J 

_ U m + l. 



t 22 J 

L U mJ 


where 


T, | = T 22 =cosh(krb) 

T 12 = Cc sinh(krb) (3-80) 

T 21 = C' sinh(krb) 

r and c,-_ are the propagation constant and characteristic impedance, respectively, b is the 
length of the element between m+1 and m. In the ZKTL code, the Zwikker and Kosten 
low frequency solution of the propagation constant for shear wave numbers below 120 
and the Kirchhoff “wide-tube” solution for shear wave numbers greater than 120 are 
used. The shear wave number is 



(3-81) 


The Zwikker and Kosten propagation constant and characteristic impedance are 
respectively 



where 


43 



n = 


(3-84) 



and y is the specific heat ratio. The Kirchoff “wide-tube “ solutions of the propagation 
constant and characteristic impedance for shear wave numbers greater than 120 are 


r= 




(3-85) 


(3-86) 

These transmission matrix elements were modified to account for the effect of the 
wave number shift caused by bias flow, which will effect cavity reactance. The elements 
in Eq. (3-80) become: 


T„ = (D _ e n ‘" b + D + e _rVb ) 


(3-87) 


T,,= 


iD + D 
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Pc£. 
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- j(e rkb -e -rVb ) 


(3-88) 


t 21 =- 


iNf 


^Ke^-e-^) 


P c c 


(3-89) 


T 22 =-j[D + e rk b +D"e _ril * b ] 


(3-90) 


where 
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k + = hS u 

1 + M c 


D + =l + (l + ir)M c 
D' =l-(l + ir)M c 

b = x m+l -x m 


(3-91) 


(3-92) 


(3-93) 


where x m+ , and x m are the x distance at locations m+1 and m, respectively. In this study 
these matrices model the cavity. Therefore, b is the length of the cavity L. 

3. 6 Error Criteria 


The method of least squares can be employed to establish the error criteria for the 
model predictions. Let E, and P, be a set of experimental and predicted results, 
respectively. Assume that the following relationship exists between Ej and P, n 

Pi = SEj + e ; (3-94) 

where e ( is the random error associated with experiment i and S is an arbitrary constant. 
Minimizing the square of the errors yields the linear regression of Ncstimatc on E where 


S = 




i = l 


(3-95) 
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and N cst , mat c=SEj. Note that this regression curve has been forced through zero. S 
measures the ratio of predicted results to experimental results. If the ratio is equal to one 
then the experimental values equal the predicted results. Therefore, the deviation of this 
ratio from one is indicative of the least squares precision error between predictions and 
experiments as seen from Fig. 3-3. 



Figure 3-3. Sketch of least squares fit between predictions and experiments. 


It is important to note that the least squares slope between predictions and 
experiments is not the same as taking the mean of the predictions divided by experiments. 
This is expressed in Eq. (3-96) 


~ F 

S * (3-96) 

N ' 

where N is the number of points in the Prediction-Experiment plane of Fig. 3-3. The 
slope in the least squares sense is more sensitive than in the arithmetic sense, and 
therefore errors are larger using least squares procedure. For example, a slope of 0.9 in 
the least squares sense does not mean that on average the predictions are within 90% of 
the experiments. On an average sense the predictions may be within say 95% of the 
experiments. 
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Another measure of error used is the mean difference between experimental and 
predicted impedance results. To calculate the mean difference, the average predicted 
impedance in the frequency range is subtracted from the experimental counterpart as 
shown in Fig. 3-4. This difference is given in terms of pc difference, since the impedance 
is normalized with respect to pc. 

Experiment 

■x . 

(Mean Impedance Difference) 

X 

Prediction 


Frequency 

Figure 3-4. Sketch of difference between prediction and experiment. 

The correlation coefficient is indicative of whether the predictions and 
experiments are following the same trend. This quantity is given by 


0 

O 

C 

03 

"O 

0 

Q 

E 


X(e,-eXp,-p) 

V i=i 


(3-97) 


where, as before, E and P stand for experimental and predicted results, respectively. The 
correlation coefficient goes from -1 to +1 if the data is negatively or positively 
correlated, respectively. Note that for the correlation coefficient the regression curve has 
not been constrained to pass through P=E=0. This was done to avoid confusion, since the 
correlation coefficient is normally presented in the literature unconstrained. 
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3 . 7 Experimental Fit to Bias Flow Impedance Model 


The impedance in the PBFIF model (Eq. (3-36)) can be written as 

k = a iG, +a 2 G 2 + ... + a n G n (3-98) 


where ai, a 2 , a n have values of one. Gi, G 2 , G n are the nondimensional group in 
that equation. For example in the PBFIF model G| and G 2 are 


16ot 

caC D d 2 


G 2 = 2.82 



caC D d 


(3-99) 


To fit the model to the experiments an optimal set of a’s is desired that minimizes the 
error between experiments and predictions. To this end an optimization procedure was 
developed and implemented. The optimization procedure was a constrained minimization 
routine that minimized an objective error function based on the slope error criterion 15 . 
Therefore, the setup for the optimization procedure was the following 


minimize — » Objective Function = abs(l - Slope) 
subject to (constraints) : 


0<a, <3 
0 < a 2 < 3 
0 < a 3 < 3 


(3-100) 


where abs stands for absolute value and Slope is the slope error criterion described in 
Section 3.6. The optimal set of a* only varied ai through av This was done since bias flow 
only affects the resistance in the models developed in this chapter. This optimization 
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routine was used for the PBFIF model and the corresponding results are shown in 
Chapter 8. The fit was performed using the 5, 10 and 15% POA samples and all flow 
rates except the highest incident bias flow rate (600 cm/s). Other fitting methods were 
tried unsuccessfully and are addressed in Appendix C. 

The ultimate goal of these impedance models is to be used to design a maximum 
absorptive liner. Appendix D addresses this ultimate goal. 
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4. Experiments 


Several perforate liner samples were tested using the NASA Langley Normal 
Impedance Tube (NIT) for the purposes of evaluating the acoustic impedance models 
discussed in Chapters 2 and 3. These samples were also tested in a raylometer for the 
purpose of determining the discharge coefficient associated with the perforate. Acoustic 
impedance was measured for the samples with and without bias flow for a frequency 
range of 1000 to 3000 Hz and sound pressure level (SPL) of 120 to 140 dB. After the 
completion of the tests, repeatability studies were performed to assess the validity of the 
experimental database. 

The normal incidence tests were performed utilizing the Non-Switching Two- 
Microphone Method (NS-TMM) 1 . This method was preferred due to the higher speed of 
data collection it offered over the Switching Two-microphone Method (S-TMM). 
Although the NS-TMM method is an established method for measuring acoustic 
impedance, the NIT facility utilized the S-TMM method rather than the NS-TMM 
method to determine acoustic impedance. Consequently, the NS-TMM method had to be 
implemented. 

In order to add bias flow to the liners, several experimental setup changes were 
required. These included the addition of a plenum chamber, the replacing of the hard 
backing plate with a high resistive mesh material, and the addition of a muffler at the end 
of the bias flow tube configuration. These and other issues involved in the testing of 
perforates in the presence of bias flow will be discussed in this chapter. 


4.1 Experimental Setup 

The NASA Langley Normal Incidence Tube (NIT) was used to make impedance 
measurements of lumped-element single-degree-of-freedom liners with bias flow (see 
Figures 4-1 and 4-2). Six acoustic drivers generate an acoustic plane-wave pressure field 
which, upon reflection from the perforate sample, sets up a standing wave along the axis 
of the 5.08-cm square tube. The perforate facesheet sample is placed at the end of the 
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tube and backed with a short 5.08-cm square cavity. This cavity is terminated with a high 
resistance fibermetal sheet designed to allow mean flow to pass through while reflecting 
almost all the acoustic signal. Three microphones are used in the test procedure. The 
microphone nearest the specimen is stationary, and is used to measure the sound pressure 
level near the surface of the specimen. Two other microphones measure the frequency 
dependent transfer functions (acoustic pressure magnitude and phase differences) 
between their respective locations. This information is equivalent to determining the 
standing wave pattern in the tube. Since the acoustic wave patterns are related to the 
surface impedance of the perforate-cavity system (test specimen), this impedance can 

2 3 

then be determined. ’ 
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Figure 4-1. Sketch description of Normal Impedance Tube with bias flow, one degree of 
freedom liner installed. 
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Figure 4-2. Picture of Normal Incident Tube. 


The surface impedance of the specimen is given by 


P 1 + R 0 . 

c = = 0 + iy 

1-R K 


( 4 - 1 ) 


where R is the complex reflection coefficient, 


R-t 

P: 


( 4 - 2 ) 


and 0 and % are the normalized resistance and reactance, respectively. 

A signal generator is used to generate discrete frequency signals that are input to 
the power amplifiers. The amplified signals are then input to the acoustic drivers. Signals 
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from the microphones are sampled and averaged using an FFT analyzer and the data is 
stored on the computer. 

Positive bias flow (blowing) is introduced through the 2.54-cm diameter inlet tee, 
shown in Figs. 4-3 and, 4-4 into a 5.08-cm square plenum chamber. The flow then 
continues through a high resistance (at least 10 pc) fibermetal sheet into the cavity 
section and through the perforate sample. The flow is exhausted through the muffler 
depicted in Figs. 4-1 and 4-2. A reference sample was tested in the NIT before and after 
the muffler was installed. Results showed the muffler had no effect on the measured 
impedance. 

To adequately measure and control the bias flow velocity in each section, four 
pressure ports were installed along the sides of each duct section before and after each 
major pressure drop in the bias flow liner. To measure the velocity through each section, 
mass continuity and the ideal gas equation are used: 


Pk^k^k — Pk+l^k+l^k+l 


Pk 


A. 

RT 


(4-3) 


Here the index k indicates the section number (see Figure 2) and p , V , P , R, and T 
are the density, velocity, absolute pressure, ideal gas constant, and temperature, 
respectively. The mass flow is measured with a laminar flow meter upstream of the inlet 
tee. Using Eq. (4-3), the velocity is calculated in each section. 
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._ BIAS FLOW INLET 



Figure 4-3. Sketch Description of bias flow one degree of freedom liner. 



Figure 4-4. Picture of bias flow double degree of freedom liner. 
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Contamination of the intrinsic perforate impedance by a “shunt impedance” due 
to plate vibration is a recurrent problem in measurements of this type 4 . For this study, 
shunt impedance effects were clearly evident for some of the perforate samples. 
Consequently, special precautions were taken to inhibit this contamination. The measured 
impedance is always a combination of the plate mechanical impedance and the liner 
acoustic impedance, which can be modeled as parallel, lumped impedances. For most 
cases, the plate mechanical impedance is high enough, relative to the liner acoustic 
impedance, to cause minimal contamination. Near the plate mechanical resonance, 
however, it becomes a significant factor; i.e., in the range of the perforate impedance. 

Near the resonant frequency of the plate, the impedance is transitioning from a 
stiffness-dominant to a mass-dominant system. Therefore, to counter resonance behavior 
exhibited in the acoustic impedance measurements, the effective plate stiffness was 
increased. This was achieved by the addition of a post support mechanism. Figure 4-5 
shows a comparison of the acoustic impedance spectra for a single perforate sample when 
it is mounted with or without the post support mechanism. The unsupported plate (no 
post support) spectra shows the resonance frequency behavior, with a drop in impedance 
above the resonance frequency. The addition of the post support eliminates the structural 
resonances of the plate (perforated sheet sample) below 3kHz; thus, the resultant spectra 
is uncontaminated by plate resonance behavior in the current frequency range of interest. 

In typical aircraft applications, this desired stiffness is achieved by permanently 
bonding a cellular honeycomb to the perforate sample. However, the honeycomb walls 
and the bonding agent cause perforate hole blockage. Since the purpose of this study was 
to study the effects of bias flow on the perforate, this blockage was unacceptable. Also, 
keeping the perforated plates unbonded allowed simple interchange of test materials. In 
contrast, the post blocked no more than one hole and it accounted for only 1 .25% of the 
total cross-sectional area of the cavity. 

Figures 4-3 and 4-4 show the post support mechanism. The post was centered 
through the fibermetal termination face into the cavity until it pushed against the 
perforate sample. A thin nut was installed on one side of the fibermetal to secure the 
fibermetal firmly. 
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S a m pie 62 


Sound Pressure Level = 120 dB 



Figure 4-5. Example of structural resonance affecting liner impedance. 

4.2 Repeatability and Error Estimation 

To ensure data quality, Perforates 52-54 and 62 were used to conduct repeatability 
tests. These samples were each tested four times over the frequency range of 1300 to 
2200 Hz (100 Hz increments) for SPL’s of 100, 120 and 140 dB. These four perforates 
were used to represent the repeatability error for the full range of percent open areas 
being tested, with no bias flow. Thus, while providing helpful information regarding the 
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NS-TMM method, these results do not offer proof of the quality of the results acquired 
with bias flow. 

For an individual sample at a fixed SPL, the mean acoustic resistance (similarly, 
for acoustic reactance) was computed from four measurements at each frequency. The 
percentage deviations (PD’s) from the mean were then computed for each of the four 
measurements. Thus, for ten frequencies at four measurements per frequency, this gave 
40 PD’s. These 40 values of PD’s were used to compute a global standard deviation for 
the selected sample and SPL, using 


X(x;-x ) 2 


s = 



N-l 


(4-4) 


where jc, is the individual PD and N=40. Since the focus of this analysis was to quantify 
the repeatability (random) error, x was set to zero; i.e., the systematic error was ignored. 
The total error from the mean (% of data lying within 95% of the mean), which is ±2s, is 
provided in Table 4-1 for each sample and SPL. 

A comparison of the Gaussian probability distribution with the measurement data 
(PD’s discussed above) is shown in Fig. 4-6. 5 This figure shows the percentage of data 
lying below a certain mean for both the ideal Gaussian distribution and the measured 
data. Clearly, the distribution is “near” Gaussian in nature; thus, computing the 
repeatability error using the Gaussian mean and standard deviation should be sufficient 
for characterization of this data. Figure 4-6 also shows that 95% of the data is within 
±7.25% of the mean. 

Table 4-1 provides repeatability data for all of the samples, at each of the three 
SPL’s tested. All twenty-four data sets show similar evidence of “near” Gaussian 
distributions of data. The measured data are shown to be off the mean value by a 
maximum of 7%. It should be noted that only 32 averages are sampled by the FFT 
analyzer for each microphone signal. In order to reduce data uncertainty, the number of 
averages could be increased. Regardless, with the data given in Table 4-1, it is reasonable 
to assume that overall measurement error is at most ±7%. 
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Figure 4-6. Comparison of ideal Gaussian distribution and data distribution (SPL of 
lOOdB, 5 POA). 


Table 4-1. Data acquisition repeatability percent error 
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4.3 Sample Description 

Each perforate was specially fabricated for this set of experiments. The samples 
are 6.35-cm square perforated plates, with rounded comers to conform to the NIT sample 
holder. The perforate’s orifices were created by the punched method. The geometric 
parameters (plate thickness, hole diameter, and percent open area), as depicted in Figure 
4, were varied for these plates over the respective ranges of 

0.24 mm < d < 1 .48 mm 
0.51 mm < t < 1.02 mm 
0.9% < POA < 16.5% 


Table 4-2 gives the target and measured dimensions for the perforate liner 
samples. The target dimensions were chosen such that only one perforate dimension was 
varied at a time. Due to fabrication inconsistencies, the measured dimensions are slightly 
different from the desired values. The numbers quoted in the table represent an average 
of several measurements, with standard deviations being within ±2% of the means for 
each measured set. The ranges of perforate dimensions were chosen to encompass what is 
typically seen in current aircraft engine liners. Several groups of 5-15 POA perforates 
with constant plate thickness and hole diameter were selected. One group of 1-5 POA 
perforates was also selected. 

A microscope was used to measure the individual hole diameters. Thirty holes 
were examined for the initial perforate. An analysis of those results concluded that only 
ten holes needed to be measured for successive samples. Perforate sheet thickness was 
measured using a micrometer, and the POA was determined by multiplying the number 
of holes in the perforate by the cross-sectional area per hole, then dividing by the total 
area of the sheet (5.08-cm square) exposed to the acoustic field in the NIT. As shown in 
the table, the fabrication process was better for the 5-15 POA perforates than for the 1-5 
POA perforates. 
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pQA = nunloer of holes x Argg of hole 
Total Exposed Area of Sanple 


Figure 4-7. Geometric parameters of liner samples. 


Table 4-2. Target and measured dimensions for fabricated perforated sheets. 


Sample 

Number 

Hole Diameter 

Sheet Thickness 

POA 


Target 

(mm) 

Measured 

(mm) 

■ Igigj 

m 

Target 

(mm) 

Measured 

(mm) 

Target 

Measured 

43 

0.024 

0.610 

0.644 

0.025 

0.635 

0.635 

5% 

5.7% 

44 

0.024 

0.610 

0.625 

0.025 

0.635 

0.660 

10% 

10.5% 

45 

0.024 

0.610 

0.599 



0.686 

15% 

14.5% 

46 

0.035 

0.889 

0.922 

0.025 

0.635 

0.635 



47 

0.035 

0.889 

0.917 

0.025 



10% 

10.7% 

48 

0.035 

0.889 

0.909 



0.610 

15% 

15.6% 

49 

0.040 

1.016 

1.039 

0.025 

0.635 

0.635 

5% 

5.3% 

50 

0.040 

1.016 

1.030 

0.025 



10% 

10.2% 

51 

0.040 

1.016 

1.029 

0.025 

0.635 


15% 

15.5% 

52 

0.024 

0.610 

0.658 



0.991 

5% 

5.9% 

53 

0.024 

0.610 

0.637 



0.965 

10% 

10.9% 

54 

0.024 

0.610 

0.639 

0.040 

1.016 

1.016 

15% 

16.5% 

55 

0.040 

1.016 

1.052 

0.040 



5% 

5.4% 

56 

0.040 

1.016 

1.037 

0.040 

1.016 

1.016 

10% 

10.3% 

57 

0.040 

1.016 

1.047 



0.965 

15% 

16.0% 


60 






























































































58 

0.055 

1.397 

1.481 

0.040 

mu 


5% 

5.7% 

59 

0.055 

1.397 

1.466 

0.040 

1.016 

1.016 

10% 

11.0% 

60 

0.055 

1.397 

1.464 

0.040 

1.016 

1.016 

15% 

16.6% 

61 

0.010 

0.254 

0.240 

0.018 




0.9% 

62 

0.010 

0.254 

0.262 

0.018 

0.457 

0.508 

2% 

2.2% 

63 

0.010 

0.254 

0.262 

0.018 

0.457 

0.508 

3% 

3.1% 

64 

0.010 

0.2541 

0.291 

0.018 

0.457 


■ 

5.2% 

65 

0.010 

0.254 

0.281 

0.018 



5% 

6.1% 


4.4 Non-Switching Two-Microphone Method 


The NIT facility has traditionally utilized a switching two-microphone method 6 
(S-TMM) that involves acquiring transfer function data between two microphone 
locations. The transfer functions between the two microphones are measured before and 
after the microphone positions are very accurately swapped by the usage of a rotating 
microphone plug. Appropriate averaging of the two readings eliminates the effects of any 
magnitude and phase differences between the two microphones. When this method is 
used with a discrete frequency source, the microphones must be swapped for each source 
frequency. While this eliminates the need for a separate calibration procedure, it is 
inefficient for the testing of a large number of test specimens. Thus, a modified version of 
the standard Two-Microphone Method 6 ’ 7 , which does not require microphone switching 
during the test process, was used to significantly shorten the acquisition duration. This 
technique will be referred to as the Non-Switching Two-Microphone Method (NS- 
TMM). 
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Figure 4-8. Magnitude and phase calibration constants for the non-switching method 
measured over a period several days to one month apart. 

Proper implementation of the NS-TMM method requires accurate amplitude and 
phase calibration, for each microphone, across the entire frequency range of interest. To 
accomplish this, the plug containing two measurement microphones was rotated such that 
the microphones were positioned in a plane perpendicular to the duct axis. For 
frequencies below cut-on for the first higher order mode, measured amplitude or phase 
differences between the microphones are due to inherent differences between the 
microphones and signal conditioning. To account for these differences, the averaging 
process of the S-TMM method was used to acquire calibration constants at each 
frequency. These calibration constants were then used in the NS-TMM impedance 
determination method. 
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Figure 4-8 depicts the variability of the calibrations over an extended period of 
time. The magnitude calibration constants vary little from day to day, but over the course 
of the test there was a variability of approximately 0.1 dB. The phase calibration 
constants have somewhat more variability (approximately 0.3 degrees). Thus, for 
improved quality, calibration constants were acquired daily for the experimental database 
presented here. 

Figure 4-9 shows a comparison of results acquired with the NS-TMM and S- 
TMM methods for a typical perforate liner. The results are almost identical. In fact, 
repeatability tests to be discussed later demonstrate more variability than shown here. 
Thus, the NS-TMM was determined to be acceptable for the current tests. 
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4.5 High Resistance Fibermetal 

For the purposes of this study, it was important to design the experiment such that 
the bias flow effect on perforate samples could be analyzed with locally-reacting acoustic 
liner models. To achieve this, one of the key elements of the bias flow liner is the 
termination at the back of the cavity. For passive liners, the termination face for the 
acoustic wave in the cavity is a highly reflective surface. To add bias flow, this 
termination must be permeable while maintaining high reflectivity. A high resistance 
fibermetal was chosen to achieve this condition. Fibermetal is a dense mesh of metallic 
strands pressed and bonded together. 

Premo also used this approach, specifically applying a backing sheet with a 
nominal flow resistance of 190 cgs Rayls at 105 cm/s. There was no mention of whether 
this value was verified; however, his results showed 0.3pc resistance difference from the 
hardwall measurement, which indicated the backing layer resistance, was not large 
enough to adequately simulate a hardwall. 

Four different methods were used to evaluate the high resistance fibermetal sheet 
used in this experiment. In most acoustic liner models, the flow resistance (sometimes 
referred to as the direct current, or DC, flow resistance) is assumed to be equal to the 
acoustic resistance at low frequencies. As described below, the first three methods 
determine the acoustic resistance using complex acoustic pressure measurements in the 
normal incidence impedance tube. Of these, the first is an indirect method, which requires 
that the acoustic resistance be educed from measurements of multiple configurations. The 
other two methods allow the acoustic resistance to be determined directly. The last 
method uses a raylometer to measure the flow resistance. Appendix E contains a 
description of the methods. 

The measured flow resistance of the fibermetal using the raylometer was 1200 cgs 
Rayls at 105 cm/s. This is much larger than the manufacturer’s quoted value of 550 cgs 
Rayls. It is also substantially larger than 190 cgs Rayls, which was the flow resistance of 
the material used by Premo in a similar test. 

The acoustic impedance measurements of the fibermetal consistently showed the 
acoustic resistance to be around 20 pc. Thus, it was expected to provide sufficiently high 
acoustic reflection to simulate a rigid termination. Figure 4-10 shows the acoustic 
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resistance measured for a selected sample liner (Sample 43) with the high resistance 
fibermetal termination versus that measured with a hardwall (highly reflecting) 
termination. The high resistance fibermetal termination causes the acoustic resistance to 
be slightly higher in magnitude (~0.05 pc) than that with the hardwall termination. This is 
significantly lower than the 0.3 pc error observed by Premo. The measured acoustic 
reactance (not shown) of the sample liner was unchanged for each type of termination. 


Sample 43 

<POA ~ 5%, t = 0.025" (0.0635mm), d = 0.024" (0.061mm)> 
SPL ■ 100 dB 



Hz 

[-♦—High Resistance Fibermetal Termination HardwaB Termination | 

Figure 4-10. Measured difference between the hardwall termination and the high 
resistance fibermetal termination. 

4.6 Cd determination 

One parameter that plays an important role in all of the prediction models is the 
discharge coefficient, Cd- The discharge coefficient is defined as the product of the 
coefficient of contraction and the coefficient of velocity. The coefficient of contraction is 
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the ratio of the area of the vena contracta to the orifice area. The coefficient of velocity is 
the ratio of the ideal to the actual velocity of the vena contracta. The vena contracta is the 
minimum flow area of a jet formed by contraction of the streamlines at the point where 
streamlines become parallel. 

The discharge coefficient can be determined by measurements of the flow 
resistance of a perforate sample compared with the flow velocity through it. The main 
device for measuring flow resistance vs. velocity is a raylometer (See Fig. 4-11). The 
pressure change caused by the sample is measured and plotted vs. the measured flow rate. 





Figure 4-11. Picture of FITF Raylometer. 
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To determine the best discharge coefficient that fits the raylometer data a least 
squares error minimization can be applied. This least squares minimization is between the 
resistance measurements and the particular model. Let y, and v, be a set of resistance and 
velocity values measured by the raylometer denoted by the subscript i. Furthermore, let ft 
be our predicted resistance model at the same point. The GE resistance model (Eq. (2-1)) 
without grazing flow in cgs units becomes 


p 3.2ixt 

d >. — V : t~t + — 7 — 

1 2000c 2 Cp ad C D 


(4-5) 


where |l, p, a, and C D , are the viscosity, density, percent open area, and discharge 
coefficient, respectfully. Equation (4-5) can be rewritten in the following form 


<K = 



(4-6) 


where ci and C 2 are defined as 


3.2p.t 
1 “ ad 2 


2000a 2 


(4-7) 


The square of the errors in resistance between the GE model and each measured 
resistance e, is 


e, =(y,-4>,) 2 

The sum of all these errors is given by 


(4-8) 
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In order to find the discharge coefficient that minimizes these errors, the derivative with 
respect to Co in Eq. (4-9) must be set to zero which produces 
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Rearranging Eq. (4-20) yields 
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Solving for the roots of this equation and making sure the second derivative of Eq. (4-10) 
is positive, produces the optimal value of the discharge coefficient that will satisfy the 
GE model for a set of experimental raylometer results. 

The GE model in Eq. (4-5) is a good approximation if the density within the holes 
does not change very much. A more refined approximation takes into account the density 
changes due to pressure and temperature. For dry air assuming an ideal gas the density is 
related to the pressure and temperature through the following relationship 8 

p = 24.0213| (4-12) 

where p, b, and T are the density, pressure and absolute temperature respectively. The 
units of p, b, and T are kilograms per cubic meter, psi, and Kelvins respectively. If Pascal 
instead of psi units for the pressure is used, then Eq. (4-12) becomes 
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p = 0.003484- 

H T 


(4-13) 


Inserting Eq. (4-12) into Eq. (4-5) produces 


24.02 13b, 3.2[it 

^ “ Vi 2000Ta 2 Cp + ad 2 C D 


(4-14) 


where the fluid (air) is assumed to have constant temperature. Equation (4-14) can be 
rewritten as 


<t>, 



(4-15) 


where ci and C 2 now are given by 
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(4-16) 


Assuming the error e s to be Eq. (4-8) and following the same optimization procedure 
outlined in Eqs. (4-9) through (4-11) yields the root equation 
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Solving for the roots of this equation and making sure the second derivative of the 
equivalent of Eq. (4-10) is positive, produces the optimal value of the discharge 
coefficient that will satisfy the GE model for a set of experimental raylometer results 
assuming varying density. 
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5. No Bias Flow Models Impedance Results 

Figure 5-1 shows the experimental vs. the PBFIF (GE-Rohr) model impedance 
results for a representative sample with no flow. When there is no flow, the GE-Rohr 
model and the PBFIF model are the same model. The numerical results approach the 
experimental results as the sound pressure level is increased in both the resistance and 
reactance. The mean slope and correlation error criteria for both the resistance and 
reactance show continuous improvements as the sound pressure level is increased (see 
Figs. 5-2 and 5-3). This indicates that the nonlinear term in the PBFIF model, which is 
associated with the acoustic particle velocity, better models the experiments than the 
linear component of this model. 


t"0.99mm, poa-5.9%> 




Figure 5-1. Experimental impedance vs. PBFIF (GE-Rohr) model numerical predictions 
for a sample with POA=5.9. 
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Cmtetidn 


Resistance Error Criteria vs. 3PL for POA«3% 




Figure 5-2. PBFIF model resistance error criteria for POA=5. 




Reactance Error Criteria vs. SPl for POA««% 



SPL 


Figure 5-3. PBFIF model reactance error criteria for POA-5. 

For the resistance (see Fig. 5-2), the mean slope error criteria starts at about 0.8 
and increases to nearly one with increasing SPL. The mean correlation starts at nearly 
zero at 120 dB and continually increases up to about 0.85 at 140 dB. The spread of the 
slope error criteria with respect to the various t/d’s increases from nearly zero at 120 dB 
to within 15% at 140 dB. For the reactance (see Fig. 5-3), the mean slope error criteria 
stays relatively constant at about 0.7 for all SPLs and the spread with respect to t/d is 
about 5%. The correlation is nearly equal to one for all SPLs with no spread in the data 
with respect to t/d. 
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«J«0.66mm, W).99mm, POA-5 9%> 




Figure 5-4. Experimental vs. GE model numerical impedance predictions for a sample 
with POA=5.92. 

Figure 5-4 shows the experimental vs. the GE model impedance results for a 
representative sample with no flow. Both the GE and PBFIF model numerical results 
approach the experimental results as the sound pressure level is increased in both the 
resistance and reactance. For the resistance (see Fig. 5-5), the mean slope and correlation 
error criteria for both the resistance and reactance show continuous improvements as the 
sound pressure level is increased. The slope error criterion starts at about 0.6 for 120 dB 
and increases to about 0.8 at 140 dB. Comparing this result to that of the PBFIF model, 
indicates that the PBFIF model is about 20% better at predicting the experimental results 
with a mean slope error criteria of about 1 for the PBFIF model vs. 0.8 for the GE model 
at the higher SPL’s. 
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Figure 5-5. GE model resistance error criteria for POA-5. 

The reactance for the GE model numerical impedance predictions (see Fig. 5-6) 
fares better than the PBFIF model with a slope error criteria of nearly 1 vs. 0.7 for all 
SPLs. It is important to note that for both the GE and PBFIF models the reactance 
prediction is the sum of the perforate and cavity reactance. Consequently, since the cavity 
reactance is much larger in magnitude than that of the perforate, the reactance prediction 
is in large part due to the ability of the transmission matrices in ZKTL to model the 
cavity. 
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Figure 5-6. GE model reactance error criteria for POA=5. 

The next couple of figures are impedance contour plots of t/d vs. frequency at 
various percent open areas. It is important to note that each t/d point in the t/d-frequency 
plane is a sample. Each of these samples had specified geometrical parameters such as 
percent open area. Nevertheless, each of these geometrical parameters varied because the 
specified and actual geometrical parameters were not the same. Besides percent open area 
the, Cd also varies from sample to sample because a single C D is calculated for each 
sample and these values are given in Appendix G. Therefore it is important to recognize 
that as t/d increases other geometrical parameters also vary, because of manufacturing 
tolerances. 

Another word of caution, realize that t/d is the part of the second nondimensional 
group in the PBFIF model (see Section 3.1). These contour plots are provided to evaluate 
how the experimental and model impedance trends against this nondimensional group. It 
is important to realize that these groups in the model are not independent of each other. 
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For example, in the PBFIF model varying t/d affect groups one, two, four, and five (see 
Section 3.1). For more information on this topic see Appendix B. 

Figure 5-7 is the impedance contour plot of t/d vs. frequency at 10% percent open 
area with no bias flow for the experiment. This figure shows how changes in t/d and 
sound pressure level, respectively affect the impedance. At low sound pressure levels 
changes in the resistance seem to be dominated more by changes in t/d rather than 
frequency, except at the high end of the frequency range. As the sound pressure is 
increased the opposite effect occurs. The resistance is now more affected by changes in 
frequency than t/d. 

These effects were expected. At low sound pressure levels the impedance is 
dominated by the linear component of the resistance, which is a function of t/d. As the 
sound pressure level rises, the impedance becomes dominated by the nonlinear 
component of the resistance. Since the acoustic particle velocity is a function of 
frequency, the resistance becomes more dependent on frequency changes. 

This figure also shows two distinctive regions of high resistance, at around 1400 
and 3000 Hz, as sound pressure level is increased. The reactance for all three sound 
pressure levels is increasing as a function of frequency, yet it is relatively constant as a 
function of t/d. This effect is due to the cavity reactance, whose resonance and anti- 
resonance occur at 2000 and 4000 Hz, respectively. 

Figure 5-8 is the GE model impedance contour plot of t/d vs. frequency at 10% 
percent open area with no bias flow. This model tends to predict similar iso-resistance 
shapes when compared to the experimental results for the lower frequency range. It also 
under-predicts the magnitude of the resistance for the entire frequency range, but the 
problem seems to be worse at the higher end of the frequency range. As the sound 
pressure level is increased, this trend seems to stay fairly constant. 

For middle and higher frequency ranges, the impedance is not predicted well. The 
model simply predicts a low resistance value through the entire region independent of 
frequency or t/d. The reactance, in contrast, is predicted fairly well for the entire 
frequency range. 
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Figure 5-7. Experimental impedance contour plot. 
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Figure 5-8. GE impedance contour plot. 
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These results were expected, since the resistance in the GE model was derived 
assuming DC or Poiseuille flow (zero frequency) through the perforate. Therefore, it is 
expected that the GE model would approximate the experimental results well in the lower 
frequency range. The higher frequency effects in the perforate’s holes, which are ignored 
in the GE model, are lost, and consequently produce a low predicted resistance value 
throughout the middle and higher spectral ranges. 

Figure 5-9 is the PBFIF (GE-Rohr) model impedance contour plot of t/d vs. 
frequency at 10% percent open area with no bias flow. The resistance predicted in this 
model approximates fairly well the experimental results for the low and middle 
frequencies for this range of t/d. This model tends to predict better the magnitude of the 
resistance especially at the higher frequency range than the GE model. 

The reactance in the PBFIF model is in better agreement than the GE model. The 
PBFIF and experimental reactances are higher than that predicted by the GE model. This 
result was anticipated, as discussed in section 2.1.2, due to the effect of viscosity in the 
mass reactance. 
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Figure 5-9. PBFIF (GE-Rohr) model impedance contour plot. 
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6. Incompressible Bias Flow Model Impedance Results 

M h is the effective Mach number (incompressible flow) at the perforate hole. It is 
defined as 


M H = V,h t 2Vbh (6-1) 

2c 

This Mach number equation is extracted from the PBFIF model’s nonlinear term (Eq. (3- 
36)). The problem with using this equation is that the actual v a h is unknown as far as the 
impedance models are concerned. In fact that is what the impedance models are trying to 
find. The experimental v a h could be acquired from the experimental results and an 
assumed relationship between the incident and perforate hole acoustic particle velocities. 
Although this method would be the most accurate, a simplified approach suffices for this 
study. 

It is of interest determine what level of bias flow is required such that changes in 
Mh due to SPL changes do not affect the overall acoustic impedance. Consequently, for a 
plane wave, the relationship between the acoustic particle velocity and pressure can be 
used to determine an approximate magnitude of v ah . Therefore, the relationship between 
the acoustic particle velocity and the SPL is 




( 6 - 2 ) 


Using the continuity principle for an incompressible fluid to relate the incident velocity 
and the hole velocity yields 


Nh| 


1 
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2caC D 

pc 


+ 2v„ 


(6-3) 
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where |M h | is used to indicate the magnitude of the Mach number in the perforate. 

The error criteria used for the resistance, was the slope and difference criterion, 
respectively. The correlation criterion was not used for the resistance, since resistance 
results showed little variation with respect to frequency. Remember that the correlation 
measures the percent change (variability) in the experiments that can be explained by the 
numerical results. Consequently, due to the small scale of variations, the correlation 
results will be of little use. The reactance, on the other hand, does have significant 
variation with respect to frequency, and therefore the correlation is meaningful. 

The discharge coefficient used for the model predictions was measured using a 
raylometer and the values are provided in Appendix G using the procedure outlined in 
Section 4.6. A discussion on the appropriateness of this procedure for determining C D 
and a sensitivity analysis is provided in Appendix F. Sample prediction results of the 
PBFIF model is presented in this chapter. For a complete database of including all the 
sample results see Appendix H. 

Figure 6-1 shows a sample experimental and PBFIF model impedance at 130 dB 
and 5.9% POA for various bias flow velocities. The figure shows increasing resistance 
with increasing bias flow rates. The resistance seems to stay constant with frequency 
except for the high bias flow rates, where it tends to decrease with increasing frequency. 
The reactance increases with frequency due to the cavity reactance. The experimental 
reactance tends to decrease as bias flow is increased. The model was not able to account 
for this effect. 

Figure 6-2 shows the slope and difference error criteria for the resistance for 
various t/d samples. The flow rates are given in cm/s and Mh (Hole Mach Number) for 
the NIT and hole velocities, respectively, on the x-axis. The slope error criterion shows a 
mean slope starting at one for no flow then decreasing to 0.7 and increasing back to 0.9 
as flow rate is increased. Therefore, the resistance numerical and experimental results are 
staying within 30% of each other in a least squares sense. 

The mean difference starts at nearly zero and increases to about one and decreases 
slightly after that. The maximum difference is about 2.5pc occurring at around 300 cm/s 
for t/d=0.689. These errors are significant; the model’s performance leaves a lot to be 
desired for this percent open area case. 
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The model underpredicts the resistances for all bias flow rates. At the higher flow 
rates the slope error improves because the experimental resista; decreases with 
increasing frequency at these flow rates (see Fig. 6-1). Since the model underpredicts the 
resistance, the decrease with frequency places the experimental resistance closer to the 
model results, and therefore decreases the measured slope error. 


sp-eztd-o.esmm, t-o.wmm, va *1.6, poa»6.9%> 
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Figure 6-1. Experimental impedance and PBFIF model numerical predictions for a 
sample with POA=5.9. 


Figure 6-3 shows the slope and correlation error criteria for the reactance of the 
same t/d samples. The mean reactance slope starts at 0.7, goes to 1.1, and then decreases 
to negative values as bias flow is increased. The mean correlation on the other hand stays 
close to one at the lower bias flow rates and then sharply decreases to 0.65. 

Note that there is significant scatter in Figs. 6-2 and 6-3 for the various t/d 
samples. This scatter is low at no flow, but significant for nearly all the flow rates. This 
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result indicates that there is some sort of interaction between t/d and bias flow velocities 
not accounted for by the PBFIF model. 

Jing and Sun 1 provide possible explanations for the scatter in error associated with 
various t/ds, and the impedance values at higher bias flow rates. They developed a model 
that accounted for bias flow. 


Resistance Error vs. Bias now Rate tor POA«5% 



Figure 6-2. PBFIF model resistance error criteria for POA-5. 


Jing and Sun 1 developed a boundary element bias flow model from Euler’s 
Equation, where the Kutta condition was applied to generate unsteady vortices due to the 
bias flow. Their model showed that different t/d and hole Mach numbers had a significant 
impact in the resistance and reactance. Their model successfully predicted a decrease in 
reactance as bias flow was increased. Since the PBFIF model does not take into account 
vortex creation in the perforate’s holes, it fails to predict the decrease in reactance, as 
well as some of the interaction between hole Mach number and changes in t/d. This 
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possibly explains the scatter seen in the error plots with t/d, which suggests that the 
models are not accounting for some physical phenomena associated with t/d. 


Reactance Error criteria vs. Bias Flow Rate tor P0A«5% 



Bias Flow Rata 


Figure 6-3. PBFLF model reactance error criteria for POA=5. 

The PBFIF model improves its prediction as POA is increased. Figure 6-4 shows 
a sample experimental and PBFIF model impedance at 130 dB and 10.9% POA for 
various bias flow velocities. The figure shows the experimental and numerical results are 
closer at 10% than at 5% (Fig. 6-1) POA. 

The error criteria for both the resistance and reactance show significant 
improvements as seen in Figs. 6-5 and 6-6. For the resistance, the mean slope starts at 
about 0.92 at zero flow, increases to about 1.18, then decreases to about 0.82. There is a 
general tendency of improving results as t/d is increased. This result was expected, since 
increasing t/d makes the modeling of perforate holes as tubes more appropriate. The 
difference error criteria for the resistance starts at nearly zero for no flow, and increases 
steadily to about 0.5pc at 600cm/s bias flow rate. Comparing the 10% vs. the 5% POA 
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difference error plots; there is significant reduction in scatter due to varying t/d for the 
10% POA case. 

SP-53<3“0.64mm, t»0.97irm, tW "1.5, POA"10.9%> 



Frequency (Ht) 


Figure 6-4. Experimental impedance and PBFIF model numerical predictions for a 
sample with POA=10.9 

The reactance error criteria shows the mean slope starting at around 0.8 at zero 
bias flow rate and increasing steadily to around 1 as the bias flow rate is increased. The 
scatter around this mean with respect to t/d increases as bias flow is increased. 
Nevertheless, the same general trend exists. That is, the slope error for the resistance 
improves as t/d is increased. The reactance correlation also shows some slight 
improvement; especially at the 600cm/s bias flow rate. This result is probably due to the 
lower bias flow rate in the perforate hole at 10% vs. 5% POA. 
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Figure 6-5. PBFIF model resistance error criteria for POA=10. 

There is a possible explanation for the PBFIF model’s improved performance as 
POA is increased. As POA is increased, the effective bias flow Mach number in the 
perforate hole decreases. For example, the highest Mach numbers achieved with the 1 5% 
and 5% POA samples were 0.126 (see Figs. 6-8 and 6-9) and 0.357 (See Figs. 6-2 and 6- 
3), respectively. The PBFIF model performs better at lower flow velocities, because the 
formation of unsteady vortices is a function of the fluid velocity, given constant viscosity 
and hole diameter. At low flow velocities, vortices are not formed, and the perforate’s 
hole physics follow the principles outlined in the PBFIF model. As the bias flow velocity 
in the hole increases, unsteady vortices start forming and the assumptions and principles 
outlined in the PBFIF model start to break down. At higher perforate bias flow rates a 
model such as that outlined in Ref. 1 should be used. 
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Reactance Error Criteria vs, Bias Flow Rate for POA*MO% 
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Figure 6-6. PBFIF model reactance error criteria for POA-10. 

Figure 6-7 shows a sample experimental and PBFIF model impedance at 130dB 
and 16.5% POA for various bias flow velocities. Comparing the sample resistance 
impedance plots for the 10% (Fig. 6-4) and 15% (Fig. 6-7) POA samples, the results look 
very similar. The resistance prediction seems to do no better in approximating the 
experimental results. The reactance, on the other hand, does show some improvement 
when going from 10% to 15% POA. 
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Experimental [Exp.] and PBFIF Model Numerical (Hum.) Results at 130 dB 



Figure 6-7. Experimental impedance and PBFIF model numerical predictions for a 
sample with POA=16.5%. 


The error criteria plot for the 15% POA samples confirms the visual trend 
comparison between the experimental and predicted results. Figure 6-8 shows the 
resistance error criteria for various t/d samples and bias flow rates. The mean slope 
criteria starts at around 0.9 at zero flow, increases to 1.2 as bias flow is increased, and 
decreases to 0.8 as bias flow is increased further. On average the numerical and 
experimental results stay within 20% of each other in a least squares sense. This trend is 
very similar to the 10% POA samples resistance slope error criteria. The resistance 
difference error criteria for the 15% POA samples starts at nearly zero at no bias flow and 
increases to about 0.3pc at 600cm/s bias flow rate. This result is slightly lower than the 
10% POA resistance error difference. The spread in difference error with respect to 
different t/d is also larger at 10% than at 15% POA. 
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The reactance slope criteria for the 15% POA samples starts around 0.8 and 
steadily increases to one, indicating that the reactance curves are within 20% of the 
experiments in a least squares sense. The reactance correlation is about one for all flow 
rates indicating a nearly perfect correlation between the experimental and numerical 
results. 


Resistance Error vs, Bias Flow Rate for POA«i5% 
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Figure 6-8. PBFIF model resistance error criteria for POA=15. 
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Reactor** Error Criteria vs, Bias Flow Rate for POA-l 5% 
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Figure 6-9. PBFIF model reactance error criteria for POA=15. 

The next couple of figures are impedance contour plots of t/d vs. frequency at 
various percent open areas. It is important to note that the each t/d point in the t/d- 
frequency plane is a sample. Each of these samples had specified geometrical parameters 
such as percent open area. Nevertheless, each of these geometrical parameters varied 
because the specified and actual geometrical parameters were not the same. Besides 
percent open area the Cd also varies from sample to sample because a single Cd is 
calculated for each sample and these values are given in Appendix G. Therefore it is 
important to recognize that as t/d increases other geometrical parameters also vary, 
because of manufacturing tolerances. 

Another word of caution, realize that t/d is the part of the second nondimensional 
group in the PBFIF model (see Section 3.1). These contour plots are provided to evaluate 
how the experimental and model impedance trends against this nondimensional group. It 
is important to realize that these groups in the model are not independent of each other. 
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For example, in the PBFIF model varying t/d affect groups one, two, four, and five (see 
Section 3.1). For more information on this topic see Appendix B. 

Figure 6-10 shows the experimental impedance for varying t/d at 100 cm/s bias 
flow for 15 percent open area. This figure shows increasing spans of constant resistance 
throughout the frequency range, especially at the lower sound pressure levels. Also note 
that increasing the sound pressure level from 120 to 130 dB does not change the iso- 
resistance distribution in the t/d and frequency plane. Therefore, bias flow has linearized 
the resistance with respect to sound pressure level at these dB levels. 

These results should be expected. Looking at the Perforate Bias Flow (PBF) 
model in section 3.1, we note that the addition of bias flow increases the nonlinear 
component of resistance. Therefore, increasing the bias flow component makes the 
nonlinear component more significant than the linear component, which depends on t/d. 
Furthermore, bias flow makes the nonlinear resistance term more dependent on the bias 
flow component than the acoustic particle velocity component. This produces the dual 
effect of reducing the resistance frequency variation, since the bias flow rate is constant 
through the frequency range, and making the resistance independent of sound pressure 
level, because the bias flow velocity is larger in magnitude than the acoustic particle 
velocity. 

Figure 6-11 shows the experimental impedance for varying t/d at 100 cm/s bias 
flow for 5 percent open area. Comparing Figs. 6-10 and 6-11 the resistance becomes 
more constant throughout the t/d and frequency range as percent open area is decreased. 
Furthermore, the resistance is now independent of sound pressure level for all three 
levels. 

The reactance does not seem to change very much between zero and 100 cm/s 
bias flow. Nevertheless, the reactance seems to be affected more by changes in percent 
open area. The 5 % open area reactance experiences a jigsaw shape with increases in t/d, 
while the 15 % open area reactance is constant with respect to t/d. 
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Figure 6-10. Experimental impedance contour plot. 
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Another trend to note is the effect of decreasing the percent open area at a 
constant flow rate. As the percent open area is decreased, bias flow becomes more 
significant. Decreasing the POA increases the effect of the nonlinear component of 
resistance in both the experiment and prediction. This is due to the l/o factor in the 
nonlinear component vs. the l/o factor in the linear component in all of the 
incompressible models. 

Figure 6-12 shows the PBFIF model impedance for varying t/d at 100 cm/s bias 
flow for 15 percent open area. Comparing this model to the experimental results, the 
PBFIF model tends to over-predict the magnitude of the resistance. Furthermore, the iso- 
resistance shapes somewhat follows a similar trend to the experimental results, yet 
improvements are needed. The effect of sound pressure level changes on this model is in 
better agreement with the experimental results. This model tends to overstate slightly the 
effect of sound pressure level changes on the resistance. The reactance is in better 
agreement with the experimental results. 
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Figure 6-12. Perforate Intermediate Frequency model impedance contour plot. 
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The PBFIF model includes a bias flow correction given in Section 3.1. This 
correction, which was derived using a frequency domain approach, is similar to that 
developed by Premo using a time domain approach (See Section 2.2.2). Figures 6-13 and 
6-14 show resistance comparisons between the models for 5% and 15% POA sample 
perforates. The reactance is not shown because bias flow only affects the resistance in 
both Premo’s and the PBFIF model. These figures show almost identical resistance 
values for all bias flow rates. 
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Figure 6-13. PBFIF vs. Premo’s impedance model predictions at POA=5.9. 
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Figure 6-14. PBFIF vs. Premo’s impedance model predictions at POA-16.5. 


1. Jing, X. and Sun, X., “Effect of Plate Thickness on Impedance of Perforated Plates 
with Bias Flow,” A1AA-99-1877. 
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7. Compressible Bias Flow Model Impedance Results 

In Section 3.2, a compressible impedance model was derived and called the 
PCBFIF model. This model was developed for compressibility effects as the hole Mach 
number approached one. Figs 7-1 and 7-2 show resistance comparisons between the 
PBFIF and PCBFIF model predictions for samples of 5.9 and 16.5 POA. Only the 
resistance is plotted, because the reactance model is the same for both models. The 
symbol with the solid line is the PCBFIF model while the symbol alone is the PBFIF 
model. These plots are plotted for increasing hole Mach numbers. The hole Mach number 
was calculated using Eq. (3-44) derived in Section 3.3 rather than Eq. (6-1), due to the 
fact that Eq. (3-44) accounts for fluid compressibility while Eq. (6-1) does not. 
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Figure 7-1. PBFIF vs. PCBFIF impedance model predictions at POA=5.9. 
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Figure 7-2. PBFIF vs. PCBFIF impedance model predictions at POA=16.5. 

The resistance results indicate that the PBFIF and PCBFIF models predict nearly 
the same results for the low hole Mach numbers. As the hole Mach number increases, the 
resistance results start diverging, with the PCBFIF model predicting higher resistance 
results. This trend is independent of the sample POA. As the hole Mach number 
increases, fluid compressibility effects become more prevalent. Fluid compressibility can 
occur due to pressure variations and entropy variations. For the PCBFIF model an 
isentropic relationship was used, therefore only pressure variations can cause 
compressibility in this model. For incompressibility the following relationship must hold 

(_V.v) = -^ = 0 (7-1) 

v p Dt 

Using the momentum equation and ignoring entropy variations, unsteady effects, viscous, 
and body forces produces 
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(7-2) 


(- V • v ) = = -M 2 \ v • V-U 2 

p Dt c 2 2 

As M becomes small, the left side of Eq. (7-2) also becomes small, and approaches the 
incompressibility relation in Eq. (7-1). Conversely, as M becomes large, Eq. (7-1) no 
longer holds and the fluid behaves as a compressible fluid. 



Figure 7-3. Incompressible and Compressible Bias Flow Hole Mach Numbers vs. 
Incident Bias Flow Rates at 5%, 10%, and 15% POA using Cd=0.76. 


Figure 7-3 shows the calculated bias flow hole Mach Number vs. incident bias 
flow rate using compressible and incompressible flow assumptions for 5%, 10%, and 
1 5% open areas, respectively. The incompressible hole Mach Number is calculated using 
Eq. (6-3), ignoring the acoustic term. The compressible hole Mach Number was 
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calculated by solving Eq. (3-44) through an iterative method. The discharge coefficient 

(C D ) used for these calculations was 0.76. 

This figure shows that the predicted hole Mach numbers using compressible and 
incompressible flow assumptions are the same for the low incident bias flow rates for 
three percent open areas as expected. As the incident bias flow rate increases, the 
incompressible and compressible flow assumptions start diverging in terms of terms of 
their prediction of the hole Mach number. This divergence starts occurring at around a 
hole Mach Number of 0.35 for all three percent open areas. From this figure a hole Mach 
number of 0.35 corresponds to an incident bias flow rate of 450, 875, and 1250 cm/s for a 
5, 10, and 15% open area perforate, respectively. Therefore, except for the highest bias 
flow rate and the lowest percent open area, fluid compressibility is not an effect. 
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Figure 7-4. Incompressible and Compressible Bias Flow Hole Mach Numbers vs 
Incident Bias Flow Rates at 1% through 5% POA using C D =0.76. 
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Samples 61 through 65 are common septum perforated plates for multi-degree-of- 
freedom liners. These septum plates range from 1% to 5% POA. Figure 7-4 shows the 
calculated bias flow hole Mach Number vs. incident bias flow rate using compressible 
and incompressible flow assumptions for 1% to 5% open area samples. As seen from this 
figure, fluid compressibility is significant at these POAs. The test range of these samples 
had to be limited to low bias flow rates to avoid choked conditions in the perforate’s 
orifices. Notice that as the POA is decreased, not only does fluid compressibility occur at 
lower flow rates, but also the rate of divergence between the incompressible and 
compressible hole Mach number predictions is greater. Therefore analysis of these 
samples must be done using the compressible bias flow model (PBCFIF) rather than the 
incompressible model (PBFIF). 
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Figure 7-5. Experimental impedance and PBFIF model numerical predictions for a 
sample with POA=0.9. 
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For the analysis of these samples the discharge coefficient used was 0.76 rather 
than the calculated C D given in Appendix G. The reason for this is that the fit used to 
determine Cd (see Section 4.6) assumes incompressible flow and therefore does not apply 
for these samples. More discussion on Cd may be found in Appendix F. 

Figure 7-5 shows a sample experimental and PBFIF model impedance at 130 dB 
and 0.9% POA for various bias flow velocities. This figure shows a rapid increase in 
resistance with increasing bias flow. The experimental resistance also shows a decreasing 
trend with frequency at the higher bias flow rates. The prediction’s agreement with the 
experimental results deteriorates as the bias flow is increased. The resistance model 
results follow the same trend as the experiments of increasing resistance with increasing 
bias flow, but tend to over-predict its values over the frequency range. 
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Figure 7-6. Experimental impedance and PBFIF model numerical predictions for a 
sample with POA=2.2. 


105 



The experimental reactance results show a decreasing trend as bias flow is 
increased. This trend is more pronounced as the frequency is increased. The model 
predictions do not account for this effect. This effect was seen before (although less 
accentuated) for higher percent open areas (5%), and more discussion on this is provided 
in Chapter 6. 

Figures 7-6, 7-7, and 7-8 show sample experimental and PBFIF model impedance 
results at 130 dB and 2.2%, 3.1%, and 5.2% POA, respectively, for various bias flow 
velocities. As with the 0.9% POA sample, the resistance results for these samples show 
deteriorating agreement between model predictions and experiments as bias flow is 
increased. The 2.2% and 3.1% samples show better agreement between the predictions 
and experiments than the 0.9% and 5.2% samples. 
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Figure 7-7. Experimental impedance and PBFIF model numerical predictions for a 
sample with POA=3. 1. 
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These higher percent open area samples also show the trend of decreasing 
experimental resistance with frequency at the highest bias flow rates. This trend was also 
seen in Chapter 6 for the 5% open area sample. It is important to note that this effect 
occurs usually only at the highest flow rate and seems to be independent of the resulting 
resistance value. Therefore, this problem may be due to the loss of backpressure in the 
experiments during the test at these low POA and high bias flow rates and not an inherit 
physical phenomena. Consequently, the model predictions will not account for this effect. 
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Figure 7-8. Experimental impedance and PBFEF model numerical predictions for a 
sample with POA=5.2 
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8. Experimentally Fitted Bias Flow Model Impedance Results 

In Section 3.6 a method of experimentally fitting (correcting) the numerical model 
to better approximate the experimental results was described. In this chapter, the PBFIF 
model is fitted using the method described in that section, with the slope error criteria 
used as the optimizing objective function. Only the resistance terms in the PBFIF model 
were fitted to the experiments. The resulting fitted model is 




caC D d' 


coC D d 


+ i 


kt _ __ V am t 8 
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aC r 


caC D d 


3k \|/(a) 


This model will be called the Fitted Perforate Bias Flow Intermediate Frequency 
(FPBFIF) model. The similarity between the PBFIF and FPBFIF models is striking. The 
optimal multiplication constant to fit the model to the experiments were 1, 0.4, and 1.2 
for the linear dc (first term), linear ac (second term), and nonlinear (third term) resistance 
terms, respectively, as seen in Eq. (8-1). Most of the correction predicted by the 
optimization routine was in the second term, reducing its influence in the model by nearly 
60%. 

Figure 8-1 shows a sample experimental and FPBFIF model resistance at 130 dB 
and 5.9% POA for various bias flow velocities. Only the resistance is plotted since the 
reactance of the FPBFIF and PBFIF model is the same. The reader is referred to Chapter 
6 for experimental and predicted reactance plots. The FPBFIF model shows good 
agreement between experimental and predicted results at the low bias flow rates. This 
agreement, as with the PBFIF model, deteriorates as the bias flow rate is increased. 
Nevertheless, the FPBFIF model seems to be closer to experimental results than the 
PBFIF model results shown in Chapter 6 and Appendix H. 

Figure 8-2 shows the error in terms of slope and difference between the 
experimental and predicted results for all nominal 5% open area samples using the 
FPBFIF model. The error plots indicate that the worst slope criterion is about 0.65, but 
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the mean hovers between 0.9 and 1 excluding the highest flow rate (600 cm/s). This is a 
significant improvement when compared to the PBF1F model, where the worst point and 
mean slope is about 0.5 and 0.7, respectively (see Chapter 6 or Appendix H). The spread 
of the slope with respect to different t/d is significant, although reduced when compared 
to the PBFIF model. There is a about a 15% improvement in terms of worst point slope 
between the FPBFIF and PBFIF model. 
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Figure 8 - 1 . Experimental and FPBFIF model impedance for a sample with POA-5.9. 

The maximum mean difference between the experiments and the FPBFIF model 
is about 2pc, although the mean is generally within 0.5pc. Comparing these results to the 
worst difference in the PBFIF model (about 2.5pc), the FPBFIF model performs about 
20% better. Although there is a significant improvement in model agreement with 
experimental results using the FPBFIF model, there is still considerable improvement 
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required. The proper modeling of the combination of low percent open area with bias 
flow is still a problem that even with an experimentally fitted model remains aloof. 
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Figure 8-2. FPBFIF model resistance error criteria for POA=5. 

Figure 8-3 shows a sample experimental and FPBFIF model resistance at 130 dB 
and 10.88% POA for various bias flow velocities. The results in this plot show significant 
improvement between the experimental and numerical results. As with the 5% open area 
sample, the agreement between numerical and experimental results deteriorates as bias 
flow increases. Nevertheless, the deterioration rate is significantly reduced for the 
nominal 10% vs. 5% open area sample. 
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Figure 8-3. Experimental and FPBFIF model impedance for a sample with POA=10.9. 

Figure 8-4 shows the error in terms of slope and difference between the 
experimental and predicted results vs. bias flow rate for nominal 10% open area samples 
with various t/d’s using the FPBFIF model. The mean slope starts at about 0.9 for zero 
bias flow, increases to 1.3 at 100 cm/s, then decreases to about 1.0 as the bias flow rate is 
increased to 600 cm/s. The spread about the mean is about 0.4 (or +/- 0.2 about the mean) 
and is fairly constant for all the flow rates. In general, as the t/d is increased the slope 
approaches one. This is expected since the perforate holes were modeled as tubes. 
Compared to the PBFIF model (see Chapter 6) the FPBFIF model fares slightly worse, 
especially for the 100 cm/s bias flow rate. 

Figure 8-4 also shows the difference between experiments and predictions. The 
mean difference tends to increase with increasing bias flow rate. It starts at nearly zero 
for no flow to about 0.5pc for the highest bias flow rate. The spread in the difference 
(maximum minus minimum difference) with varying t/d also follows this trend of 
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increasing with bias flow increase. The spread starts at nearly zero and increases to about 
lpc (or +/-0.5pc about the mean). 


Resistance Error vs. Bias Flow Rata for POA»io% 
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Figure 8-4. FPBFIF model resistance error criteria for POA=10. 

Figure 8-5 shows a sample experimental and FPBFIF model resistance at 130 dB 
and 16.48% POA for various bias flow velocities. The results at the for the 15% open 
area sample are similar to the 10%, and both show better agreement than the 5% open 
area sample. 

The slope error criteria of the 15% nominal POA is very similar to that of the 10% 
POA (see Figs. 8-4 and 8-6), although the spread about this mean seems to be reduced for 
the 15% vs. 10% POA. The mean difference error criteria for the nominal 15% POA is 
less than 0.1 pc for all bias flow rates, and less than 0.05pc excluding the highest bias 
flow rate (600cm/s). The spread of the difference error is also within 0.1 pc (+/- 0.05pc) 
excluding the highest bias flow rate (600cm/s). 
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Comparing these results with the PBFIF model (Chapter 6), both the slope and 
difference error criteria for the nominal 15% POA samples are very similar. The FPFIF 
model performs somewhat better than the PBFIF model, especially at the higher bias flow 
rates. For example, the mean slope error at 600cm/s bias flow rate is about 0.8 for the 
PBFIF vs. 0.9 for the FBFIF model, a 10% improvement in the least squares sense. 

SP-*4<d««,84mm, tM "16, POA«1«.5%> 



Frequency (Hrj 


Figure 8-5. Experimental and FPBFIF model impedance for a sample with POA- 16.5. 
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Bta Flow R»to 

Figure 8-6. FPBFIF model resistance error criteria for POA=15. 
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9. Conclusions and Future Work 


There were four objectives to this study: 

1. Acquire a quality experimental database of educed normal incidence 
impedances for perforates plates with and without bias flow 

2. Develop improved impedance models for perforated plates in the presence of 
bias flow 

3. Evaluate the models developed against the experimental database 

4. Use the experimental database to improve the model’s performance and 
evaluate the resultant model performance 

To achieve objective one, normal incident impedances were experimentally 
educed from perforated samples in the presence of bias flow. These samples were tested 
in the NASA Langley Normal Incident Tube (NIT), where several experimental issues 
were addressed to test these samples. 

These issues included the modifications to the NIT to introduce bias flow, the 
repeatability of the experiments, the implementation of a non-switching two-microphone 
method, mechanisms to eliminate the perforates plate resonances, the accurate 
determination of bias flow speeds through the perforates, and the impedance 
determination of the high resistance fibermetal between the plenum chamber and the 
cavity. Figures 4-1 and 4-3 show the NIT experimental setup. Having addressed these 
issues, a quality database of educed perforated plate impedances was acquired, and 
therefore objective one was achieved. 

Objective two was addressed by first determining the previous impedance 
modeling methods. These methods mainly involved the determination of the impedance 
of perforated plates without bias flow. Three impedance models, which included the 
effect of bias flow, were developed. Two of the impedance models were developed by 
modeling the perforate orifices as infinite tubes with the inclusion of end corrections. 
These models assumed incompressible and compressible flow assumptions, respectively, 
between the far field and the perforate orifice. The incompressible and compressible 
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models were called the PBFIF (Eq. (3-36)) and PCBF1F (Eq. (3-52)) model, respectively. 
These models were implemented and their results were evaluated against experiments. 
The third model was developed following a similar approach to Hersh99. This model was 
not implemented and more information on this model may be found in Appendix A. 

The PBFIF and PCBFIF models were evaluated against experiments to address 
objective three. The frequency range for which this model was evaluated was from 1000 
to 3000 Hz. The PBFIF model performed better for higher (15%) rather than lower (5%) 
percent open area (POA) samples. Based on the slope error criteria, the model’s 
resistance predictions were on average within 20% and 30% for the higher and lower 
POA’s, respectively. The mean difference between the experiments and predictions 
started at nearly zero for the low flow rates and increased to about lpc and 0.5pc for 5% 
and 15% open area respectively, at the highest incident bias flow rate of 600cm/s. The 
numerical model predictions of lower POA perforates contained more scatter, with 
respect to varying t/d’s, than at higher POA. This was true regardless of whether the 
slope, difference or correlation error criterion was used. The scatter also tended to 
increase as the flow velocity was increased for any given POA. 

Results at zero flow indicated the linear terms required more correction than the 
nonlinear terms for the PBFIF model. When this model was compared to GE’s model, 
which lacked the linear resistance frequency dependent term, the PBFIF model 
outperformed the GE model when compared to experiments in the middle and higher 
frequency range. 

As flow was added the nonlinear term became more important than the linear 
terms. The resistance of the perforate became more independent of t/d, SPL, and 
frequency as bias flow was increased (perforate linearization). The point where this 
linearization occurred was a function of the bias flow rate and percent open area. For 
example, the 15% POA samples experienced linearization at lOOcm/s incident bias flow 
velocity for the SPL’s of 120dB and 130dB, but not 140dB, while the 5% POA sample 
was linearized for these SPL’s at this incident bias flow rate. 

The experimental reactance results showed a decreasing trend as bias flow was 
increased. This trend was more pronounced at the lower percent open areas (<5%) and 
higher frequency range. The PBFIF model was not able to account for this effect on the 
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reactance. A possible explanation for this shortcoming in the model was that this model 
did not account for the formation of unsteady vortices. Another model by Jing and Sung 1 
that accounted for this effect was able to predict the decrease in reactance as bias flow 
was increased. The PBFIF model showed good agreements with Premo’s bias flow 
correction. 

The PCBFIF model and the PBFIF predicted nearly the same results for low hole 
Mach numbers. As the hole Mach number was increased the PCBFIF model predicted 
higher resistance values than the PBFIF model. It was shown that the predicted hole 
Mach number was the same using the compressible and incompressible flow assumptions 
except for the highest bias flow rate (600cm/s) and lower percent open area (5%). 

The septum perforates with POAs between 1% and 5% were also analyzed. The 
results showed that the rate of divergence between the predicted hole Mach number using 
compressible and incompressible flow assumptions increased as the POA was lowered. 
Therefore, the PCBFIF model was used to analyze these samples. The septum perforate 
results showed decreasing agreement between model results and experiments as bias flow 
was increased. The 2.2% and 3.1% samples show better agreement between the 
predictions and experiments than the 0.9% and 5.2% samples. 

To address objective four, a fit of the PBFIF model was performed using the 
experimental data. The fit was performed using an optimization routine that found the 
optimal set of multiplication constants to the non-dimensional groups that minimized the 
least squares slope error between predictions and experiments. The fitted model indicated 
that most of the correction required to better approximate the measured results was in the 
frequency dependent linear term of the resistance. The fitted model results showed 
significant improvement over the original model for the low percent open area 5% 
samples. The fitted model performed about the same as the original model for the 10% 
POA sample and slightly better for the 15% POA sample. 

There are several potential research areas for future work in the modeling of bias 
flow acoustic perforated plate liners. The determination of a suitable grazing flow model 
that accounted for bias flow was only preliminarily addressed in this dissertation. There 
are several other parameters that have future research potential in the area of bias flow 
acoustic liners. These include the effects of broadband acoustic signals and temperature 
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variation, to name a few. Another important area of future research is the determination 
of Cd- Appendix F discusses some of the issues associated with the methods used to 
determine Cd in this study, and provides a sensitivity study of the effect on the resistance 
results due to the uncertainty in Cd- There are also several experimental fluid mechanical 
visualization studies that could be performed to complement the development of other 
theoretical models. 

1. Jing, X. and Sun, X., “Effect of Plate Thickness on Impedance of Perforated Plates 
with Bias Flow,” AIAA-99-1877. 
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Appendix A: Perforate Boundary Layer Bias Flow (PBLBF) 
Model 

The Perforate Boundary Layer Bias Flow model is based on the principles 
outlined by the Hersh99 model 1 . The original model includes cavity reactance and needs 
to be modified for this study since ZKTL calculates the cavity reactance in a separate 
module. Hersh99 utilizes a control volume approach, where they write the conservation 
of vertical momentum within the holes. 



Perforate Hole 
v eh> S s 


Perforate Hole Wall 


Figure A-l. Control volume for perforate hole. 

The momentum equation from Fig. A-l for the modified model is 

p hS) H^-+p h [27CrvLdr-pv^S 1 = P,S, -P,(S, -S s )-k w S w (A-l) 
dt J 0 

where R is the radius of the perforate orifice. Hersh uses a similar equation to Eq. (A-l) 
except that he includes the cavity reactance. He then groups terms, simplifies, and 
rewrites that equation in terms of semi-empirical parameters, which produces 

p h S,H^^- + p h S 1 f^--l jv| +K vis v e = P s S s (A-2) 
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where the parameters G, H, and K v i s are determined semi-empirically. Hersh 
demonstrates that G can be related to the viscous boundary layer profile by 


G = 


^BL (1 K BL ) 

0-s.J s, nv 


where Kbl is given by 


^ _ V BL^BL 

*^BL 


(A-3) 


(A-4) 


Let v c be a harmonic acoustic disturbance with a superimposed mean bias flow 
velocity of the following form 


v e =v a e i0) ' +v b (A-5) 

and let 

v e = l v a + 2v b v 0 Je'“* + v|; (A-6) 

where Eq. (A-6) forces nonlinear and linear oscillations to be of the same frequency. 
Let’s further assume that 

p s =Pa e '“ I+ P b (A-7) 

Substituting Eqs. (A-5) through (A-7) into Eq. (A-2) yields 

ip,S,Ho>v.e“ +P.S0- 1 )v.e“ + v„ 'f + K * (v,e“ + v b )= (p,e“ + p.Js, (A-8) 

Rearranging and collecting terms produces 
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[(iA + B + 2v b )v a + v a 2 - Cp a Js‘“‘ + [v 2 + Bv b - Cp b Je° = 0 (A-9) 

V ^ ' ' 


=0 


"Y" 


where 

(A- 10) 
(A-ll) 

(A- 12) 

In order for Equation (A-9) to be satisfied the coefficients of the terms in that equation 
must be zero. For the first coefficient this produces 

(iA + B + 2v b )v a + v a -Cp a =0 (A* 13) 

Multiplying Eq. (A- 13) by v a produces 

iAv 2 + Bv a + 2v b v a + v a = Cp a v a (A-14) 



Integrating both sides over the period yields 



(A- 15) 


121 


where v] has been replaced by |v a |vf in the fourth term to ensure that it remains 

positive. Assuming a harmonic solution for v a of the form Cos(cox) and simplifying 
yields 


iAv L + Bv L + 2v b v^ s + — 2V2v^ = ^Cpcv; 

371 


,2 

rms 


(A- 16) 


Dividing both sides by C pcvj^. produces 


S = — p t B + 1 1 v m» + 2v b + iA] 
pcC 


(A-17) 


Variables A, B, and C depend on K V i S , H, and G. These variables must be related 
to the geometrical parameters of the perforate that produces them. Hersh developed these 
relationships in his paper through semi-empirical means. These relationships are 


K.„ = 


rcdp 




(A- 18) 


K =13 + 50.35 - 




(A- 19) 


K =3.2 + 1.95 


( t T 1 07 


V' 


(A-20) 


The parameter H becomes 


H = [a + (1 - a)e _bP “" Ji e 


(A-21) 
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where 



b = 1.537 + 6.532 - - 


(;M;) 


+ 0.179^ — 
d 


d c = t + 0.85d 


I 



The parameter G becomes 


G it -nPnon 

= a[l-e J 


where a and b are redefined as follows 



11.876 + 4.48 


b = • 


1 + 0.15 


dll 


(A-22) 

(A-23) 

(A-24) 

(A-25) 

(A-26) 

(A-27) 

(A-28) 
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These semi-empirical relationships assume a viscous boundary layer profile with an 
inviscid core in the perforate’s hole. The boundary layer thickness has been assumed to 
be equal to the diameter of the perforate hole. 


References 
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Impedance Model,” ALAA 99-1825, 1999. 
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Appendix B: Non-dimensional Analysis and Functional 
Determination 

The impedance of liners can be analyzed numerically and experimentally by using 
principles of non-dimensional analysis. Let the experimentally determined non- 
dimensional impedance % be expressed as 

% = F(G,,G 2 ,...,G n ) (B-l) 


where F is a function and G,, G 2) ..., G n are non-dimensional groups. In general these 
groups are complex, and are themselves a function of dimensional geometric and 
physical parameters. 

Theoretical impedance models, such as those developed by Crandall and Kraft, 
have claimed that the function F is separable. This means Eq. (B-l) can be rewritten as 

£=F,(G 1 )+F 2 (G 2 )+... + F n (G n ) (B-2) 


The different theoretical models assume different F’s and G’s depending on the 
assumptions and approximations of each model. In general Gj, G 2 , ••■» Gn are not 
independent of each other; i.e., they may be functions of each other. Nevertheless the 
functions Fi, F 2 , ..., F n can be approximated from the experimental results through the 
following procedure. Taking the derivative of Eq. (B-2) with respect to G| yields 

3^ _ dF,(G,) | dF 2 (G 2 )3G 2 ^ , dF n (G n )3G„ (B . 3) 

3G, dG, dG 2 3G, dG n 3G, 


Rearranging and integrating both sides produces 


'-’ir 

F,(G,)= J 


3E, dF 2 3G 2 
3G, dG 2 3G, 


dF„ 3G n 
dG„ 3G 


dG, 


(B-4) 


1 J 
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where Gn and G| f are the initial and final Gj, respectively. Consequently the functions F, 
are 


w lf 

F.(G.)= 1 


dF 2 3G 2 


g l 3G, dG 2 3G, 


f 2 (g 2 )= j 

Cji L 


3^ dF, 3G, 


3G 2 dG, 3G 2 


dF„ 3G n 
dG n 3G, 

dF„ 3G„ 
dG„ 3G, 


dG„ 


dG 


2 » 


F„(G.)= | 


at 


dF, 3G, 


dF„.| 3G„_, 


[3G n dG, 3G n dG n _, 3G n J 


dG„ 


(B-5) 


If the functions Gj are independent or weakly dependent. Equation (B-5) becomes 


F,(G,)=°|'^-dG„F,(G 1 )=°|'idG 2 F.(G,)=°fidG, (B-6) 

G. d °i G. oG, j 3G„ 


Therefore experimental impedance results should be plotted against each of these non- 
dimensional groups G|, G 2> G n . 
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Appendix C: Alternate Experimental Fit to Bias Flow 
Impedance Model 

Another method of determining the functions Fj(G.) is to assume their form and 
use a least squares fit criterion. Assume the following form for F, 

F i (G i )=a i G i +e i (C-l) 


where a, and e, are an unknown constant and random error corresponding to the non- 
dimensional group i, respectively, to be determined from experimental results. Double 
subscripted notation is now introduced, where the first and second index indicates the 
non-dimensional group number and experiment number, respectively. For example G, 2 
indicates the value of dimensionless group 1 for the second experiment. In what follows 
the subscript in the single subscripted quantities indicate parameter variations. For 
example, ^ indicates non-dimensional experimental value i and a, indicates a constant 
associated with group j. Notice, that although both indices are the first indices in these 
variables, they mean different things. Letting the square of the errors be e 2 and applying 
the relationships just mentioned produces 

e 2 = (£, -a 1 G„-a 2 G 2l -...-a n G nl ) 2 + -a,G 12 -a^^ -...-a^^) 2 ^ ^ 

+ ”- + 0»m -a iG| m - a 2 G 2m — ••• -a n G nm ) 2 

for n dimensionless groups, and m experimental values. To minimize these errors, the 
partial derivative of e 2 with respect to a; is taken and set to zero. Doing this yields 


— = 2fe,-a l G tl -a 2 G 21 -...-a n G nl )G ll +2fc 2 -a 1 G l2 -a 2 G 22 ... a^X^ 

8a, 

+ ... + 2 (^ m — a,G lm — a 2 G 2m — ... — a n G nn ,y3im — ® 


Rewriting Eq. (C-3) in summation notation results in 
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m m mm 

a .Z G ij G . J +a 2X G 2jGi J +- + a n £G nj G ij =£^G a (C-4) 

J =1 j=l j=l j=l 


This is the equation for group i. For n dimensionless groups Eq. (C-4) in matrix form 
becomes: 


m m 

£ G U G ij ZGjjG.j . 

j=l J=1 

m 

ZGnjG.j 

j=l 

a i 


m 

Z^ G ,j 

j=l 

X^2jG 2 j 

J=1 j=l 

m 

•• £ G "j G 2j 

j-i 

r 2 

► “ i 

m 

I^G, 

j=l 

m m 

S G lj G nj Z G 2j G nj ' 

. J=> j-1 

m 

lG nj G nj 

j =1 

V 


m 

Z^jGnj 

>' J 


Letting G, A, and Z be the first, second and third matrices respectively, of Eq. (C-5) 
produces 

[gKa}= {z} (C-6) 


Therefore, A can be found from 


{aHgHz} (C-7) 

where [g]~ 1 is the inverse of matrix [g]. 

Using the least squares approach proved to be difficult to implement. The 
problem stemmed from the fact that the difference between largest and lowest magnitude 
dimensionless groups in [g] varied significantly (high matrix condition). Techniques 
such as singular value decomposition (SVD) 2 were tried and found to be unsuccessful. 
Therefore another fitting technique was used. 
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Appendix D: Absorption Ct f ficient 


The idea behind modeling bias flow acoustic liners is to be able to use these 
models to generate an optimally designed absorptive liner. In what follows, it will be 
shown that absorption coefficients, reflection factors and insertion losses are related to 
the acoustic intensity. The liner was designed using an incident plane wave analysis as 
shown in Fig. D-l. 


y 



Figure D-l. Plane Wave Model in ZKTL. 


I and R in Fig. D-l stand for the incident and reflected acoustic wave, respectively. M 
stands for the grazing flow Mach number and note that the e" 10 ” convention is now being 
used for this section only. Using continuity of particle displacement at y=0 results in the 
following relation for the reflection coefficient, R f , 


R f 


£(l + M sin <(>)cos <|) - 1 
£(l + M sin <J>)cos<J> -t- 1 


(D-l) 


where 



(D-2) 
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For the normal incidence tube, design M will be zero. For each spectral component, the 
mean squared pressure is given by 



( 4jrf 




(D-3) 


where 6 is the phase of R f . Since the reflected wave is coherent with the incident wave, 
there is a mixed term as shown in Eq. (D-3) that contributes to the SPL. This causes the 
pressure field to vary spatially because of interference effects. At y=0, Eq. (D-3) 
corresponds to the input pressure for ZKTL which is the total acoustic pressure that is 
held constant during the iteration algorithm. Thus, for a given SPL, changing the liner 
parameters/configuration changes the ratio of reflected energy to incident energy but has 
no effect on the SPL values. Minimizing the reflected energy is the goal of the liner 
design procedure. 

The acoustic intensity (time-averaged energy flux) of the incident wave normal to 
the liner must be conserved which is 1 



(l + Msin<J))|l| 2 cos<|) 
2pc 


(D-4) 


The difference between the incident and reflected intensities represents the energy 
absorbed by the liner: 


The absorption 


( w »)- (UMli £ M — 6- N’) 


coefficient is defined as 


(D-5) 


= (^4 = l-|R r | ! (D-6) 

K) 
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Substituting Equation (D-l) into Equation (D-6) produces 


a _ 49(l + M sin 4>)cos<(> 

[©(l + Msin^cosifH-l] 2 + [^(l + Msin<j>)cos<|)] 2 ^ ^ 

which is the relation used in ZKTL to compute the absorption coefficient. Also, in ZKTL 
the "R.F." that is computed is equivalent to 


R.F. = l-a = |R f | 2 (D-8) 

which is the ratio of reflected energy to incident energy, (w Ry )/(\V ly ) . The insertion loss, 
"I. L.", computed by the code is 


I.L. = 



(D-9) 


Thus, an optimum liner requires maximum ot, minimum R.F. and maximum I.L. 
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Appendix E: High Resistance Fibermetal Determination 
Methods 

T ransfer 




Indirect Method 

Test 1: 

Test Sample: 

Perforate sheet 

Cavity Termination: 

Hardwall 

Bias Flow: 

Off 

Test 2: 

Test Sample: 

Perforate sheet 

Cavity Termination: 

High resistance fibermetal sheet 

Bias Flow: 

On 

Test Sample: 

Direct Method I 

High resistance fibermetal sheet 

Cavity Termination: 

Hardwall 

Bias Flow: 

Off 

Test Sample: 

Direct Method II 

Open (no sample installed) 

Cavity Termination: 

High resistance fibermetal sheet 

Bias Flow: 

On 


Figure E-l. Summary sketch of three methods used to determine acoustic resistance of 
fibermetal sheet. 
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4.5.1 Indirect Acoustic Method 

A low resistance perforate was tested in the normal incidence tube using two 
configurations. For the first configuration, the perforate had a cavity with a hardwall 
termination. For the other, the termination was the high resistance fibermetal under 
investigation, which could allow bias flow to be passed through the tube. Using the NS- 
TMM method described earlier, the acoustic impedance of the perforate sample was 
determined tor each of these configurations. 

The impedance measured at the surface of the perforated plate for the hardwall 
termination is given as £ sl . Similarly, the impedance measured with the fibermetal 

termination is £ s2 . These two impedances are the sum of the individual impedances of 
each liner element; i.e., 


5.,-5,+U 

(E-l) 

5a -5, +5* 

(E-2) 


where £ p is the perforate impedance, and £ cw and £ cf are the cavity impedances with the 

hardwall and fibermetal terminations, respectively. For this method to work £ must be 

independent of the test configuration. The only way to assure £ p is constant is for the 

perforate to be linear (independent of SPL and bias flow). Tests were conducted over the 
bias flow and SPL range of interest in this study with the high resistance fibermetal 
termination installed. The results of these tests indicated that the selected perforate 
sample was acceptably linear. 

Subtracting Eq. (E-l) from (E-2) and solving for £ cf produces 


k=U+k-5.. 


(E-3) 
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This is the cavity impedance with the fibermetal backing. The cavity impedance 
due to the hardwall is given by 


^ cw = -iCot(kL) (E-4) 

where L is the cavity length, and k the wave number. The lumped-element impedance 
of the fibermetal is related to the cavity impedance with the fibermetal termination £ cf 

by the following relationship 




where the wave numbers k : and k r are 


k: = ■ 


k. = 


1 + M 
k 

1-M 


(E-6) 


and M is the Mach number in the duct. 

The major advantage of this method is that £ sl and £ s2 can be measured 

accurately. Measurements near nulls of large standing waves are avoided by properly 
choosing the perforate material, consequently improving the accuracy of the 
measurement. Perhaps more importantly, bias flow effects on the fibermetal impedance 
can be studied. The major disadvantage of this method is that more measurements and 
calculations are required to determine the impedance of the fibermetal. 

4.5.2 Direct Acoustic Methods 

There are two ways to measure the fibermetal impedance directly. In Method 1 
the fibermetal sheet under investigation is installed as the “test sample” in the NIT, with a 
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hardwall termination. The NS-TMM method is then used to measure the normal 
incidence acoustic impedance. This method offers the advantage of requiring only a 
single measurement, and consequently is a fast method for determining the impedance 
properties of the fibermetal. Its main disadvantage is the inability to measure how bias 
flow affects the fibermetal impedance. 

In Method II the fibermetal under investigation is used as the cavity termination. 
For this method, no sample is installed; thus, the cavity is left open. The fibermetal 
impedance is determined using Eq. (E-5), where £ cf and are the measured impedance 
at the standard test plane (where sample surface is typically located) and fibermetal 
lumped-element impedance, respectively. The advantage of Method II is the single-step 
process in measuring the impedance of the fibermetal. Furthermore, this method allows 
for the determination of bias flow effects on the impedance of the fibermetal. The major 
disadvantage to this method is the potential for measuring near nulls of large standing 
waves, with the accompanying increased potential for measurement error due to large 
changes in SPL over the diameter of the measurement microphone. 

4.5.3 Raylometer Method 

Pressure 



Figure E-2. Raylometer experimental setup for determining fibermetal resistance 
(Raylometer Method). 
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The raylometer measures the DC flow resistance of the high resistance fibermetal 
backing. Figure E-2 shows the typical experimental setup utilized for this experiment. 

The non-dimensional resistance, 9, is 


e= Pl ~ Ei (E-7) 

pcV b 

where p, , p, , V b , p , and c are the pressure reading before the fibermetal, pressure 

reading after the fibermetal, velocity in the duct, fluid density, and the speed of sound, 
respectively. The flow resistance is assumed to be a linear function of velocity of the 

form 


9 = A + B— (E-8) 

c 

Several values of flow resistance versus flow velocity were acquired in this 
experiment. These values were then curve-fitted to determine A and B. 

The advantage of using the raylometer is the speed of acquiring the data. The 
disadvantage is the assumption that the resistance of a DC flow measurement is 
equivalent to the real part of the acoustic impedance. This assumption is not entirely 
correct because frequency dependence is ignored. 
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Appendix F: Dynamic C D 


Jing and Sun 1 took issue with the empirical nature of the determination of Cd- 
Their model showed that the formation of the vena contracta is a time dependent problem 
in the presence of an acoustic field, and determined an approximate steady state value of 
0.61 for C D . This value of C D is lower than the normally measured experimental value of 
0.76 or higher. Jing and Sun did not account for the presence of mean flow, and therefore 
the variation of C D with increasing bias flow is not known and is an area of potential 
future research. 

For this study, the Cd used for the model predictions, except for the compressible 
model results, was a measured raylometer Cd- This approach may be questionable but it 
has been ingrained in this field for nearly three decades. The value of C D greatly affects 
the predicted results and in general it is a function of frequency, SPL, etc. Tradition and 
the empirical nature of this field (more than physics) has generally resulted in the use of a 
constant C D in previous studies. Figures F-l to F-9 show the PBFIF model resistance 
error criteria sensitivity to C D ’s of 0.76, 0.66, and 0.86 for 5, 10, and 15 POA, 
respectively. These plots show significant variation in the results depending on the value 
of C D used. Therefore, the accurate determination of this variable and consequently an 
increased understanding of the physics of this variable in the presence of an acoustic field 
with mean flow is an important area of future research. 
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Figure F-l. PBFIF model resistance error criteria for POA=5 and C D =0.76. 
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Figure F-3. PBFIF model resistance error criteria for POA=5 and C D =0.86. 
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Figure F-4. PBFIF model resistance error criteria for POA=10 and C D =0.76. 
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Figure F-5. PBFIF model resistance error criteria for POA-IO and Cp-0.66. 
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Figure F-6. PBFIF model resistance error criteria for POA=10 and C D -0.86. 
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Figure F-7. PBFIF model resistance error criteria for POA=15 and Cd=0.76. 




Figure F-8. PBFIF model resistance error criteria for POA=15 and Cd=0.66. 
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Figure F-9. PBFIF model resistance error criteria for POA-15 and C D -0.86. 
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Appendix G: Kluge and FITF Raylometer Database 


Kluge is an assembled raylometer or “old” raylometer. This raylometer contains four 
pressure ports to measure the pressure drop across the sample. 

FITF is a high-end raylometer with 160 measurement locations, statistical measurement 
outliers elimination, and curve fitting to minimize tailing effects* near the sample. 


The pages that follow provide plots of resistance vs. velocity for the perforated 
plates samples. Each page contains three plots, where the first, second and third plots 
provide results using the Kluge, FITF, and Kluge & FITF raylometer results, 
respectively. Cda and Cdb are the calculated by doing a linear regression of the plots, and 
equating the resultant intercept and slope to (c,/C D ) and (c 2 /C D ), respectively, in Eq. (4- 
6). This creates two equations, where C DA and C D b are obtained by of solving C D in the 
resulting intercept and slope equation, respectively. Cdbst is calculated using the 
procedure outlined in Section 4.6. The Co used for the prediction was the Cdbst of the 
FITF raylometer. 


* Tailing effects is near field effects close to the sample. The pressure field starts tailing 
(curving) near the sample. 
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Resistance g*s-cm 2 ) Resistance g^s-cm 2 ) ^ Resistance g^s-cm ) 


SP-43 

<d=0.64389 mm (0.02535 In), t=0.635 mm (0.025 In), POA=5.B65%> 
Kluge Raytometer Data 
Cd =0.48302, Cd b =0.73207, Cd brt =0.7262 
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SP-44 

<<1=0.82535 mm (0.02462 in), t=0.6604 mm (a028 In), POA=1 0.473%> 
Kluge Raytometer Data 
Cd a »0. 50314, 0^=0.82184, Cd 6tt =0.78062 
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Resistance s-cm 2 ) Resistance gfls-cm 2 ) Resistance s-cm ) 


SP-45 

<d=0.59919 mm (0.02359 fci), t=0.6858 mm {0.027 In), POA=14.48B%> 
Kluge Raytomeier Data 
Cd =0.73873, C<1=0,88988, C0 M =C.B4949 



RTF Raytometer Data 
Cd g =0. 70075, 0^=0.95833, 00^=0.93929 
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Resistance I* W| ^ ResMence ,* s-cm 2 ) R.slsttnc. W eW) 


SP-46 

<d»aS2151 mm (0.03628 In), t “0.635 mm (0.025 In), POA=5.376%> 
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SP-47 


<d=0,91669 mm (0.03609 in), t=0.635 mm (0.025 In), POA=10.69%> 
Kluge Raytometer Data 
Cd s =0 45515, 0^=0,79605, 0^=0.78275 
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RTF Raytomatar Data 
Cd a =0.37361, 0^=0.88101, Cd^sO.87249 
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SP-48 

<<1=0.90832 mm (0.0358 In), 1=0.6096 mm (0.024 In), POA*15.627%> 
Kluge Raytometer Data 
Cd t *0.3761, 0(^=0.83915, 00^=0.79944 
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FITF Raytometer Data 
Cd a =041429, 00^=0.91055, Cd^aO.89129 



Comparison of Best Rt Cd for Kluge and RTF Rayfometers 
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Resistance e^SH^n 2 ) Resistance ctfe-cnr) Resistance grts-cm ) 


SP-49 

<d=1.0394 mm (0.04092 In), t=0.635 mm (9.025 In), POA=5.293%> 
Kluge Raytoroeter Data 
Cd # =-0.097045, 00^=0.61439. 0^=0.66666 
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RTF Raytometar Data 
Cd g =1.0993, Cd b =0.7577B, Cd^O.80579 



Comparison of Best Rt Cdfor Kluge and RTF Raylometers 
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Resistance gAs-cm 2 ) Resistance g/(s-cm 2 ) Resistance pAs-cm 2 ) 


SP-50 

<d=1,Q305 mm (0.04057 In), t=0.635 mm (0.025 In), POA=10.148%> 
Kluge RayVxnetw Data 
Cd 1 =a47008, Cd^sO.77935, Cd blt =0.77266 



RTF RaytomeUr Data 
Cd # =0.72281, €*,*0.7733, Cd^-0.78877 



Comparison of Best Fit Cdfdr Kluge and RTF Raytometers 
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Resistance {tf s-cm ) 


SP-51 


<d=1.0232 mm (0.04052 in), t=0.635 mm (0.025 in), POA=15.474%> 
Kluge Rayfcurteter Data 
Cd # =0.43137, Cd,=0.82203> Cd bs ,=0.8003B 



o Experimertai 
— 1st Order Fit 
-- Best Cd 






SP-52 

<d=0.65811 ram (0.02591 in), tsQ.9906 mm (0.039 in), POA=5.919%> 
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SP-53 

<430.63729 mm (0,02509 In), t=0.8652 mm (0.038 In), POA=1 0.877%> 
Kluge Raytometer Data 
Cd t =0. 50062, 0^=0,76149, Cd bgt =0.72784 
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Cd a =0.5398B, 0^=0.8286, Cd^O.81658 
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SP-54 

<d=0.63906 mm (0.02516 In), t*1.016 mm {tt04 InX POA»16.481%> 



156 








157 








SP-56 

<d=1.0373 mm (0.04084 In), t=1.016 mm {0.04 In), POA*10.283%> 
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Resistance g/fs-cm 2 ) Resistance g/fs-cm 2 ) Resistance gfls-cm 2 ) 


SP-57 

«t*1.0472 mm (0.04123 Ini, t=C,9652 mm (0038 In), POA=16.021%> 
Kluge Raytocneter Data 
Cd t =0.43778, 0(^=0.76095, Cd fcrt =0.73705 



Comparison of Best Fit Cd for Kluge and FITF Raylometers 



Velocity (cmte) 
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Resistance g fls-cnr) Resistance g^s-cm 2 ) Resistance s-cm 2 ) 


8P-68 

<d=1.4808 mm (0.05829 In), t=O.S906 mm (a039 In). POA=5.671%> 
Kluge Raylometer Data 
C^s-0.14252, Cd b *0.B7274, Cd^O.71323 



Comparison of Best Rt Cd for Kluge and RTF Raylometers 
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SP-59 

<d=1 .4656 mm (0.0577 h), !«1.016 mm (0.04 ln>, POA=10,882%> 
Kluge Raytometer Data 
Cd^.28810, Cd b =0.80709, Cd fcrt =0.7855 
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Resistance tfs-cm 2 ) Resistance tfs-cm 2 ) Resistance g*s*cm 2 ) 
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SP-60 

<d=1.4643 mm <0.05765 in). 1=1.016 mm (0.04 in), POA=1 6.641 %> 
Kluge Raytometer Data 
Cd,=0.52216, 0^=0.8337, (^,=0.62292 
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Cd a =tf.30O3$, Cd^sQ.SIHS, 0*1^=0.87076 



Comparison of Best Ftt Cd for Kluge and RTF Raylometers 
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SP-61 

«t=0 23878 mm (0.00944 In), t=0.5334 mm (0.021 In), POA=C.898%> 
Wuoe Raytom^er Data _ 


3000 


b 2000 
o> 1500 


VC 


r— -r r- 

""" ' 1 

0 



L- - 

, ! 1 1 1 

t5 

o Expertmertai 
— IstOnJerFit 
- BestCd 

i p- 

i9n 14 


20 


40 


60 



163 






400 r 


SP-62 

<d=0.26238 mm (0.01033 Ini, t=C.50B mm (0.02 In), POA=2.184%> 
Kluge Raytometer Data 
0^=0.59407, 0^=0.78635, 00,^=0.78429 
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SP-63 


<d=0.26213 mm (0.01032 in), t=0.508 mm (0.02 In), POA=3.13%> 
Kluge Raytometer Data 
Cd =0.67363, 0(^=0.69652. 0(^*0.70059 
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Resistance grt(s-cm 2 ) 


SP-64 

<d =0.291 08 mm (0.01 146 In), 1=0.508 mm (0.02 In), POA=5.242%> 
Kluge Raylometer Data 
Cd # =0.36522, CdL =0.75706, Cd =0.68328 




Velocity (cm/s) 
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Appendix H: Experimental and PBFIF Model Impedance 
Database 

The pages that follow provide plots of PBFIF model and experiment impedance 
results. These plots are not the complete database. To get the complete database contact 
Mike Jones in the Structural Acoustics Branch at NASA Langley Research Center. 
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SP-45<d»0.59919mm, t=0.6B38mm, P0A=14.4B9%> 
Experimental (Exp.) vs. PBFIF Model Numerical fNum.l Results at 0 cm/e 



Frequency (Hz) 


0 



SP-46<d“0.92131 mm, t-0.635mm, POA»5.37e%> 



Frequency (Hz) 





Reactance Resistance 


SP-47<<f*0.91 669mm, t=0.633mm, POA«10.69%> 
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Resistance 


SP-48<d=0.90932mni ) t=D.6096mm, POA=15.627%> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 0 cm/s 



Frequency (Hz) 
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SP-49<d«1 .0394mm, t=0. 633mm, POA<=5.293%> 
Experimental (Exp.) vs. PBFIF Model Numerical {Num.y Results at 0 cm/s 




Resistance 
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Reactance Resistance 


SP-51<d=1. 0292mm, t=0.635mm, POA*15.474%> 
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Resistance 


SP-C^d'sQ. 63729mm, t*0. 9652mm. POA*1 n.BT7%> 




Frequency (Hz) 
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Resistance 
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Resistance 
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SP-38<d=1 4806mm, t=0.9906mm, POA“5.671%> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 0 cm/s 
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SP-59<cM .4556mm, t»1.016mm, POA*10.9B2%> 
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SP-43<d«0.64389mm, t=0.635mm, t/d "0.96619, POA«5.665%> 
Experimental (Exp.) vs. PBFIP Model Numerical (Num.) Results at 130 dB 
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Resistance 


SP-44<d=0 62835mm, t^O.eSMmm, t/d «1.0561, POA=10473A> 
Experimental (Exp.) vs. P8FIF Model Numeri cal (Num.) Resorts at 130 dB 
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Resistance 


SP-45<d*0.59919mm, t=0.6858mm, t/d =1 .1446 POA=l4 AR9*4> 
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SP-46<d=0.92151 mm, t=0.635mm, tfd =0.68908, POA-3.376%> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 130 dB 
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Resistance 


35 


SP*47<d B 0.91 669mm, t=0.635mm, f/d *0.69271 , POA=1 0.69%> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 130 dB 



Frequency (Hz) 
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Reactance Resistance 
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Reactance 


SP-49<d=1 ,0394mm, t=0.635mm, t/d =0,61093, POA=5.293%> 
Experimental (Exp.) vs, PBFIF Model Numerical (Num.) Results at 130 dB 



2000 2200 2400 2600 2800 3000 



192 










Resistance 


SP-5CKd=1 .0305mm, t“0.635mm, tfd «0.61622 r POA=1 rtR 

Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 130 dB 
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SP-51<d=1 .0292mm, t-0.635mm, t/d -0.61690, POA=1 3.474%> 
Experimental (Exp.) vs, PBFIF Model Numerical (Num.) Results at 130 dB 
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Reactance Resistance 
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Resistant* 
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Reactance Resistance 


SP-54<d“0.63906mm, t**1.016mm, t/d *1 .5896, POA“1 6.481 %> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 130 dB 
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Resistance 


SP-58<d“1 ,4806mm, t=0.9906mm, t/d =0.66907, POA=5.671%> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 130 dB 
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Resistance 
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Reactance Resistance 


SP-6(Xd=1 4643mm, t»0.506mm, t/d *0.34692, POA=1 6.641 %> 
Experimental (Exp.) vs. PBFIF Model Numerical (Num.) Results at 130 dB 
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